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The Distribution of the t and F Statistics for a 
Class of Non-Normal Populations! 


RautepH ALLAN BRADLEY 


University of North Carolina and Virginia Agricultural Experiment 
Station of the Virginia Polytechnic Institute 
INTRODUCTION 

The study of errors of experiments involving small samples began 
with the fundamental work of “Student” (1908). A test of the hypoth- 
esis that the mean of a group of independent observations has a pre- 
specified value was derived by considering the moments of the sample 
mean and of the sample standard error on the assumption that the ob- 
servations came from a parent normal population. R. A. Fisher (1925, 
1928) extended “Student’s’” work to obtain a test of the difference be- 
tween the means of two groups of observations, and rigorously proved 
the results of “Student”, at the same time introducing the use of the 
t statistic in this test. Fisher also obtained a test of the hypothesis that 
k groups of observations, assumed to come from normal populations 
with equal variances, have equal means. The z statistic used in that 
test is a simple monotone function of the more common statistic F used 
in the analysis of variance. This work is summarized by Fisher (1935) 
in an expository paper. 

Uncertainties still exist regarding the practical validity of these tests 
when the assumptions implicit in their derivations are not satisfied. 
This matter is the subject of a sequence of papers by W. G. Cochran 
(1947), Churchill Eisenhart (1947), and M. 8S. Bartlett (1947). The as- 
sumptions in the analysis of variance are that 

(i) The experimental errors are normally distributed; 

(ii) The experimental errors have a common variance; 

(iii) The experimental errors are independent; and 

(iv) The group characteristics tested and the environmental effects 
are additive. It is with violations of assumption (i) that we shall deal 
in this paper. 

There are many published papers on the subject of the non-normality 
of experimental errors, but the available information is unsatisfactory. 
The errors introduced in the significance levels of the ¢ and F tests by 


1 Research supported in part by scholarships awarded by the National Research Council of Canada 
and by the Ontario Research Foundation, and presented to the Institute of Mathematical Statistics, 
Chicago, December 27, 1950. 
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the non-normality of parent universes require more detailed and com- 
prehensive study. Such a study is necessary not so much because these 
tests are the appropriate ones for non-normal populations (and they 
are not usually best in these cases), but because it is desirable to use 
a common analysis which will give good results for a wide class of popu- 
lations. Assurance that the results obtained are reliable is required. 

Approximate methods of checking the validity of ¢ and F tests are 
diverse. Transformations of the observed variates are discussed by 
Bartlett (1947). E. J. G. Pitman (1939) has given a general procedure 
for obtaining tests adapted to specified non-normal populations. Com- 
binatorial and rank order tests are considered by Pitman (1937), Fisher 
(1947), Harold Hotelling and Margaret Pabst (1936), and by E. 8. Pear- 
son (1931). The principal discussion of limiting distributions of the ¢ 
statistic is presented by Kai Lai Chung (1946) and his results are valid 
for large sample sizes. Verifications of the efficacy of the ¢ test are ob- 
tained by experimental sampling by Pearson (1937), W. A. Shewhart 
and F. W. Winters (1928), and T. Eden and F. Yates (1933). P. R. 
Rider (1929), Victor Perlo (1933), and J. Laderman (1939) have dis- 
cussed the exact distributions of ¢ for rather trivial special cases and 
very small sample sizes. Bartlett (1935) and R. C. Geary (1936) have 
investigated errors in the significance levels of ¢ by considering the first 
few terms of a Gram-Charlier series. Their results are approximate, 
depending on a knowledge of measures of skewness, and on neglecting 
powers of the sample size of lower order than its inverse and powers of 
the skewness of higher order than the first. 

While the theoretical results so far obtained have not been particu- 
larly useful, still mathematical analysis appears to be the method by 
which more definite confirmation of the applicability of analysis of vari- 
ance to non-normal populations may be derived. This paper was sug- 
gested by an asymptotic evaluation of corrections for non-normality at 
high significance levels by Hotelling (1947) and by an extension of his 
work by the present author (1950). 

It is generally known that the cumulative distribution functions of 
t and of F may be written in the form of multiple integrals which depend 
on as many variables of integration as the number of observations in 
the sample and whose domains of integration depend on ¢ and on F. We 
show that, for a large class of non-normal populations, the density func- 
tions of ¢ and F may be written in terms of integrals whose domains of 
integration are independent of ¢ and of F. These density functions in 
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turn are used to expand the cumulative distribution functions in powers 
of (' and of F~*. The expansions are particularly useful at high signifi- 
cance levels and it is believed that these results represent a considerable 
improvement over previous investigations of the effects of non-nor- 
mality. 

In the following sections, derivations for the one-sample ¢ test are 
given in detail. Only final results (derived by analogous methods) are 
shown for the two-sample ¢ test and the F test.” Examples are developed 
for those cases in which the observations are drawn from the Cauchy 
and the “squared hyperbolic secant” distributions. 


THE PROBABILITY DENSITY FUNCTION 


Let x, --- , tw be a sample of N independent observations from a 
population with density function, f(u). Without loss of generality, we 
use the ¢ statistic to test the hypothesis that the mean of u is zero. In 
order to obtain a suitable expression for the cumulative distribution 
function of ¢, Gy(t), it is required that f(u) have the following properties: 

(a) f(u) > Ofor all -x« <u < om; 

(b) f(u) is a continuous function of u for —«# <u < «. Now, using 
the familiar formula for t, we write 


Gy(t) - | ead Mee ies [ Tse) aes. 


. (2){v 3 — #8/(N — »}" <t 


If we use a prime to denote a column vector, we transform to new 
variables of integration by setting 


(2.2) y =Tr’, 
wherein we require that 
~ g-2 7M 0 ey 7 
67? gt? —2(6)* 


{N(N — 1)}7?{N(W — 1)}7?{N(N — 1)}7) 








(2.3) r= toe {N(N ty 1) “tat ne {(N _ 1)/N}” 
NT | ee N7M 
vee n-? nN? 


*A complete discussion is available (Bradley, 1949) in a dissertation in the library of the University 
of North Carolina at Chapel Hill. 
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Upon application of the transformation inverse to (2.2), we have 

tt NP yy — {G — 1)/i}""ya + EE + DM ys | 
+e) + ((N — DN} yy 

for all = 1,--- , N if it is understood that y is taken to be zero. 


Since T is orthogonal, the domain of integration in (2.1) is simplified 
and it follows that 
Ge = | wih, ile toneie Us I ax. 
- =l j=1 
(2.5) ” . . 
N- —1/2 
YN lems WN T< p}"" < 
We use f; to simplify the writing of f(x;) when x; is taken to be a func- 
tion of the variables of integration. 

By transforming to generalized spherical coordinates*, we may re- 
duce the domain of integration to one involving a single variable of 
integration. This procedure involves setting 

Yn = pcos 
Yi. = p sin ®, cos ®, 


Y2 = p sin ® sin & cos 


Yn-2 = P sin P; sin ®, sin ®; eas sin Py_» cos Py_} 
Yv-1 = p Sin ®, sin & sin #3 --- sin @y_» sin Py_; . 


It is easy to verify that the transformation (2.6) will be one-to-one 


when we take the ranges of integration, 0 to a for 4, ets Py», 0 to 
2x for dy_,, and 0 to « for p. The jacobian is p* sin * “@, sin 
‘~, --- sin dy_» . In addition, if we set 

(2.7) t = (N — 1)'” cot 9, 


the domain of integration in (2.5) is given by 0 < #4, < x. With the 


understanding that 9 is a function of ¢, we may now write 


Gu = [ sin *“® «| [ Lf W 5 


(2.8) 
- p* "sin *“*@, - -- sin by_s dp Tl d®; Jan. 


3 See, for example, Hotelling (1931). 
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It follows at once that 


dG y(t) a N-2 
de 


r w p2xr( N N-—1 
{| | -: | I ll ih p. ‘sin **@, --- sin by_» dp ae. | 
0 0 “2S a, 2 $\=8 


The argument of f; in (2.9) may be simplified by an inverse sequence 
of transformations. We begin by setting 


(2.9) 


(2.10) — ui = (y:/sin ®) | o,-0, fae bes Ow mf, 
The jacobian is (p** sin N~*p, .-- sin @y_») and p’ = -  “. Te 
final transformation is obtained by defining new variables, 

y= a + Ota eee + {(N — DIN) un 

eo = -7 "y+ 6 "m+: + {NN — 1)N} "uy 
(2.11) w= —2(6) "ue +s + {(N — DN} wy 


vn = —(N = 2Y(N = 1) ye + O(N — DN ee. 


This transformation, although not orthogonal, has a simple jacobian. 
Let the matrix of the transformation be written as \ and then 


1 
Mr’ = 
1 
1/N 
so that |X| = N~”. The jacobian of the transformation from the 
u-coordinates to the v-coordinates is simply N*”. Further v, + --- + 
vy = (N — 1)'?N*?uy_, , and if (2.11) be written v’ = du’, wv’ = 


ur’ru’ and 
N—1 N-—2 N-1 
2 2 y—1. 2 2 r 7-1 2 
>. v35= z uz; + N uy-1 = z uz; — (N — 1)N“uy-1. 
1 1 1 


Then, 
N—1 : N—1 a N—1 2 
(2.12) dw = Sei+ (La). 


1 
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With the help of these transformations, and upon noting that the new 
variables are unrestricted on the real axis, we have 


dG () _ _N"? gin *9 [ be fl 
co : Lae ~ 


f sin 3) \w" bs yp? + (x u) cot 8 — Leh] 


; TL s| sin on™ = ve + (XZ) cot 8 + o] 
Sat + (= So.) Maw. 


= = ng may be derived from (2.13) when we recall that 
= (N — 1)” — © from whence sin @ = {1 + @/(N — 1)}~” and 


Gy(t) = SOr(l) - 7 with d@/dt = —(N — 1)” sin *e. Hence, 





(2.13) 








Gy) = nl 7" 1 & /(N ee |) oa 


ges [slo + ¢/(N — 1)" 
{nav —1)" / = vi t+ (= ns) ~ b> ot 





(2.14) Thy| a + #/(N — 1)” 
{nv -1)" / > +(E0) + ot] 





/nN—1 as a Se 
Dit (Xo.) II ae. 
1 

We have now obtained an expression for the probability density 
function of ¢ for the large class of populations satisfying conditions (a) 
and (b). When f(u) is the normal density, it is very easy to show that 
(2.14) gives the familiar “Student” distribution. When f(u) is not the 
normal density, it may be possible to evaluate Gx(t) directly in rare 
cases and for small sample sizes. More difficult cases may be treated 
by the method of the following section. 
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A SERIFS REPRESENTATION OF THE CUMULATIVE 
DISTRIBUTION FUNCTION 

Gy/(t) will be expanded in a Taylor Series in powers of 1/t. To do this, 
it is necessary that f(w) have as high order continuous derivatives as 
terms used in the expansion. We then may require the condition, 

(ce) d* f(u)/du* is a continuous function of u for —%«% < u < © and 
& = j, 3, 3, +->. 

It will be convenient to write 


9 


N-1 N—1 2) 1/2 
(3.1) P(v) = aaa dof + (= ».)} and n=N-1 
1 1 


in the expression for dGy(t)/d9 in (2.13). In that formula, we set 
(3.2) zi = »;| cos O|, i=1,---,n, 


and this is a permissible transformation so long as 89 ¥ 2/2. The ja- 
cobian is simply | cos 8 | ~". We also define 


(3.3) @ = |tan@| = n'?|¢|~. 


When @ is between 0 and 2/2, d@/d® = cos’ 8 andt is positive; when 
@ is between 1/2 and x, d0/d@ = —cos’ 6 and ¢ is negative. Upon mul- 
tiplying (2.13) by d6/d® and applying (3.1), (3.2) and (3.3), we obtain 


dGy(t)/d® = —Ne*” ps [ baie aos a2.) | 
(3.4) > <7 


${ Pte) —@® . 2 P(2) I dz, 


for positive values of ¢, and 


dGy(t)/d® = ver [ ves [ | I f{-P@ + a2.) | 
r.) — 00 1 
(3.5) F 
J{-P@ —#) 2 P(2) I] des, 
1 1 
for ¢ negative. 

In the following discussion we shall use the Taylor Expansion of 
Gy(t) in powers of about the point @ = 0. The form (3.4) will be con- 
sidered in detail and the result for (3.5) may be written down in a 
similar way. 

dGx(t)/d® may be taken to be the product of two functions of and 
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successive derivatives may be obtained by the application of Leibnitz’ 
rule. We are concerned with the derivatives of 


(3.6) A(@) = 6", 
and 


n n . 
(3.7) B(@) = {TI f{P(2) + v2) |s{PC0 — X 2 


evaluated at 6 = 0. The required derivatives of A(#) vanish except for 


(3.8) d*~*A(®)/de*~ = Tn). 
A generalization of Leibnitz’ rule permits us to write 
; 1 ae k k! 
¢B@)/de po = = = ——_ 
(3.9) aot: + *+an,=k Gig! *** Ay! 


(—2zi)* att --- 2a" sf” {P(z)} --- fo" {P(2)}. 


By means of (3.5), (3.8) and (3.9) and, again, by the application of 
Leibnitz’ rule, we may state the result that 


N&" (n — 1 . 
Gy(t) = Gy(t) jo=0 — ns (" 0 )r@) / cee 





~ : n No" +k k ae 1 C) 

(3.10) [ B@)|»0P@ I dx—--- - = ne )r in) [ 
* PB(®) | in| |. 

Lo Ge bw PO TT des — = 


If we recall that @ = n"’ | ¢|~', when ¢ > 0, G@y(t) |e = 1, and 


: ~ = Nn"” . n—l1 —n . 
Gv) = 1 n! ro ( 0 ) . 
eee / B(®)| 9-0 P(z) I] dz; 
ne eae Nn” oo rin)(” + k ieee — [ 
(n + k)! k 0 


eae 





(3.11) 





leno P(@) TT des - ++, 


for positive values of ¢. 








If 
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When ¢ is negative, we need only define. 
(3.12) B*($) = {PO —® 2 a ITs —P(z) + &2;}, 


and use (3.5) in place of (3.4). Then, when ¢ < 0, Gy(t) | e.0 = 0, and 
the expression equivalent to (3.11) for this case is 


n/2 ta e-) 
Gd ao rin)(” ; ') | t ‘ad | 


n! 20 





— [ B*(@) leno P(z) I] dz; 
Nn? n+k— ') er 
ee | k | ¢| [. 


(3.13) 


leno P(2) I] des + ee 6 


These expansions of Gy(¢) hold promise of rapid convergence for large 
absolute values of ¢, and afford a method of obtaining at least a few 
terms of the distribution function of ¢ when the population density 
function, f(w), satisfies conditions (a), (b), and (c). We now develop a 
systematic method of evaluating the terms of these series. 


THE REDUCTION OF THE COEFFICIENTS 


The expansion of Gy(¢) depends on the integral of d*B(®)/d#* | g.oP(z) 
over unrestricted values of the real variables of integration. The k* 
derivative of B(®) may be expressed in two parts, one of which is inde- 
pendent of the density, f(u), and hence common to all functions Gy(t), 
and another much simpler expression depending on the particular popu- 
lation under consideration. The first step in such a procedure is to note 








that P(z) and >-} z; are symmetric functions of 2, --- , Zn. Hence, 
from (3.9), 
d‘ B(®) (—Dai)2ft- -- 2%" ey f"? {P(z)} 
4.1 = kin! : = 
(4.1) db le ‘: ae Mm!s**Ny! II a;! 
ayS-++S ay 


where a , -*- , @ are assumed to take on p distinct values, the first 
° ° ° 4 
of which occurs 7, times, the second n,. times and so on. 


4 We use the symbol = to read “has the same integral as’’ with respect to the variables present. Thus, 
if u = dv, J“, 9(u) du = f° g(dv)d dv, and g(u) = dg(d»). 
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Many linear transformations may be used to transform P(z) to a 
function of a sum of squares, but, in order to take advantage of the 
ordering of the exponents of z , --- , 2, , we must choose a transfor- 
mation such that the complexity of the expressions for the z’s in terms 
of the new variables decreases from z, to z, . The transformation best 
suited to this plan is one we have already used (2.2, 2.3). It is only 
necessary to replace N by n in the matrix, I’. If the new variables are 
taken as w’s, we may write 


z= nw, — 0 — 1) wii + {tlt + 1}; 


(4.2) —1/2 
+--+) + f{(n — 1)n} wr-1, a7=1,---,n. 


Now, 
( n n 2) 1/2 ( n 1/2 
2 r—1/2 2 2 r—1/2 
P(z) = {2 zia+ (> :.) >» N =<>> w; + nv’, eh 
1 . 3% ; 
. / r1/2 ° ° r—1/2 ° 
If we set w, = Nw, , we have jacobian N-”, and upon replacing 
, r1/2 Ps ° . . 
w, by w,/N° and dropping the prime, the expressions required are 
, AT\—1/2 —1/2; 1/2 for2 \)—1l/2 
2; = (nN) “uw, —7t “A — 1) "win + {it + 1)} W; 
—1/2 
tee + {Qn — In} we, 


te woe 2 
n 
/ar\l/2 
(4.3) >> Zz: = wr(n/N)”, 
1 
n 
2 r\ 1/2 
P(z) = (> wi/N)"”. 
1 
For brevity, let us write z; in terms of w, , +--+ , Wn as 2:(w). 


A transformation to generalized spherical coordinates is again speci- 
fied. A first step sets 


(4.4) Wi = pe, t=1,---,n, with 2. e, = 1. 
1 
The second step requires that we take 


€: = sin 9, sin9,_; --- sin Q 


(4.5) é; = sin 9, sin9,_; --- sin 9,4; cos 9;, 
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The jacobian of the transformation is 


(4.6) p” J(Oz, --- , On) 
with J(@., --- ,9,) = sin "0, sin "“°O,_, --- sin ®; , and the ranges 
of integration are 0 to w for 9; , --- , 9, , and 0 to 2m for ® . 

Since the transformation (4.3) is linear, z7‘ is a homogeneous function 
of w,--:, W, of degree a;, and we may write 2“) = p*i{z(e)}*. 
In terms of the new variables (4.1) may be written 

d‘B(®) | 
dt | <0 
en 
@| an n \ 

s p 2i'(e) «++ 2n"(e)\ —en ala 

(47) 2kint yO 4'O (¢) J 
5 Dr tots + an=k ——  o J (@2 oo ee é,) 
OF cxo-ae Mm! --* Np! 


prt Pp PP o/ VN) 
. j=0 a;! 
where the accumulated jacobians are included. 

It is well to note at this point that we have divided each term in the 
summation into two parts, one of which is a function of p alone, and 
therefore an even function of w,,---, w,, Whereas the other is in- 
dependent of p and is, excluding J(®., --- , 9,), a polynomial of degree 
k in terms of «,----: , €,. The previous analogous part was then a 
polynomial of degree / in w;,--- , w, , and when k is odd, at least one 
of wy, +--+ , W, appeared as an odd function in every term of 

ao 
a" (w) < —w, /u\ : 
VN) 
Since P(z) is an even function of the w’s, the integrals of d*B(®) /d@*| g_ 
P(z) are integrals of odd functions whenever k is odd and then vanish 
independent of the form of f(u). As a further observation, when k is 
even, certain terms of 


212(w) coe 2°*(w) { —W, =e 
N| 


will contain one or more of the variates w,,--- , w, with an odd ex- 
ponent and such terms may again be considered to be odd and may be 
omitted. This property carries over to 6 ,---, €, and we need only 
consider terms which consist of products of ¢,--- , €, all with even 
exponents. For example, when k = 4, we consider only terms containing 
¢; and e;¢;. Elaboration of these simplifying properties will later be in 


zi'(w) --- 2 
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order but the discussion is included here since they are more easily 
detected in terms of w, , --- , w, and carry over to &, +++ , én. 
Let us define 


; 7 rm - n! f_ TP ae 


(4.8) . 
s J (02, cre 8.) II 257 (6) 

j=1 
where a ,°-*: , @ is a combination of integers such that a + --- + 


a, = k, and the summation over the set A represents a summation over 
the fixed combination with a) taking each distinct value of the combi- 
nation once and only once with a, S --- S a, taking the remaining 
values after a) has been determined. p refers to the number of distinct 
values in this remaining set with the first occurring n,; times and so on. 
Cyle, k; a, -+* , Gn) is then well defined. We shall further write 


(4.9) Dy(0, k; ao, -** , an) = p”*™* TT (f°? (p/-/N)/a;!} 





j=0 

and note that permutations of a ,--- , dn leave Dy(p, k; ao, +--+ , Gn) 
invariant. Using (4.7), we have 

k 

3 
a @ P(z) 

dd | a-0 

(4.10) 


= Nk! p ® Cre, k; ao, ++, Gn) Dw(p, k; ao, --* 5 Gn), 


comb. 
aot+*: *+an=k 


where the summation has now been reduced to one over the combi- 
nations of integers @ + --- + a, = k instead of one over the permu- 
tations of the same set of integers. The expression for P(z), given in 
(4.3), has been included in the right hand side of (4.10). 

The integral of d*B(®)/d®*|goP(z) may be writen as the product of 
two integrals. It is only necessary to define 


(4.11) Din(k; ao, +++ 5 Qn) - | Dy(e, k; a, +++ An) dp 
and 


Cin (k; M&,°** 5 Qn 


(4.12) - © Qn 
— [ ff [ Cwle, k; do, +++, Gn) dO, --- dOe. 
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Then, 
[ eo 
wo adv 





P(z) II dz; 
o=0 1 
(4.13) 
at k! >> Cin (k; a, +++, Qn): Dink; ao, "+95 An). 


_ N comb. 
Got: + ++an=k 


In exactly the same way as that described above, it is possible to write 
™ ” d'B*(®) - 
vr pleco SA ici P " 
[. [. db $=0 (2) I * 


_kt 2 


~ = comb. Cin(k; ao, hs , On) -Diw(k; do, ee » Gn) 
N aot: + +-+onmk 





(4.14) 


where 


(4.15) Dx(o, hi; a0, +++, an) = 0” TT (f°? (—p/V/N)/a31} 
1 
and 


(4.16) Diw(k; ao, io an) = | Dx(o, k; a, gaa: » An) dp. 


D 
THE EVALUATION OF Ciy(k; ao, +++ , dn) 


Crn(k; ao, ++ , dn) need only be evaluated once for each set of values 
a + ---+ a, = k since it is independent of the particular density 
considered. As previously indicated, C;, vanishes for all odd values of k 
and we then need examine only the even values. In order to obtain 
Cry it is first necessary to write Cy as an explicit function of & ,--- , €: 
and to proceed with the transformation (4.5). We omit all terms of 
Cy which contain an odd power of any ¢;, and, furthermore, from 


symmetry, e]°! --- €2" may be replaced by ev! --- €°" where Bi pitty 
B, is a rearrangement of 6,,---, 8, with 6; <--- < 6. Thus, for 
k = 4, say, Cy may be reduced to a linear function of ¢', and exen—1 
only. 


When k = 0, Cr(e, 0; 0, --- ,0) = J(@2, +++ , On) and Cry(0; 0, ---, 
0) = i see / / sin *~*6, sin "*O,-1 --+ sin 03 dO, --+ dOs. If we 
0 0 0 


w/2 
recall that 8(m,n) = 2 sin ””' cos *” '@ d®, it is easy to show that 
0 


(5.1) Crx(0;0, +--+ ,0) = 2n"?/T(n/2) = Sar. 


This defines the constant, S,_-1. 
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The calculations of C;y are tedious but straight-forward. We illus- 
trate by evaluating Cry(4; 2, 1, 1,0, --- , 0). Cwr(e, 4; 2, 1, 1,0, --- , 0) 
is obtained by substituting for z,;(e) using (4.3), which defines z;(w), 
in (4.8). Upon applying the rules of procedure given at the beginning of 
this section, we have 


3 2 
Cy(e, 4; 2, 1, 1,0, ---,0) = J(@2, ---,0,) {es een 2 


4 2 4 
2 2 n—n — 12n+ 12 \ 
— En€n-1 ca { r a 


2nN \° 
Crw(e, 4; 2, 1, 1, 0, --- , 0) will depend on two integrals for which the 
variables of integration are @.,---, @,. It is necessary to use (4.5) 


to express ¢, and e,_; as functions of the required angles. Integration 
depends on the definition of the Beta integral and it can be shown that 


| wee | [ €),(O2, --* , Ondd (Oc, --- , On) dO, --- dOs 
0 0 ~0 
= 3S8,1/{n(n + 2)} 


and 

® © pdr 

[LL de + 0d ae, + 001@:, +++ 50.) dO, ++ ds 
0 0 “0 


= S,-1/{n(n + 2)}. 
Then, 





2 — 
(52) Crm(4;2,1,1,0,---,® — Sit (* on 2). 


~ N(n+2)\ 0 2 

Values of the C,, coefficients are listed below for values of k up to six. 
They are sufficient for obtaining the first four terms of the expansions 
of Gy(t), for populations whose densities satisfy the mild conditions 
imposed, when combined with the coefficients D;y which alone depend 
on the population density. 

We list’ the coefficients: 


(5.3) k=0 & = Sn1 
Ciw(0; 0, «+ 0) = 

k=2 & = Sar 
Cin (2; 2, 0, - + 0) = & 


§ See footnote 2. 
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Crv(2; 1, 1, 0, «=> 0) = —}6, 
k=4 & = S,4/{N(n + 2)} 

Crv(4; 4, 0, «=> 0) = 3nd, 
Crv(4; 3, 1, 0, «++ 0) = —3nd, 
Crw(4; 2, 2,0, --- 0) = 3(n® + 2)6, 
Crv(4; 2, 1, 1, 0, --- 0) = —3(n® — 3n + 2) 
Crw(4; 1, 1, 1, 1, 0, --> 0) = 3(n® — 3n + 2)d, 

k=6 ds = Sps/{N?(n + 2)(n + 4)} 
Crw(6; 6, 0, «+: 0) = 15n’é. 
Crw(6; 5, 1, 0, +++ 0) = —15n7b, 
C1w(6; 4, 2, 0, «+= 0) = 3n(n? + 4)d. 
Crv(6; 4, 1, 1,0, «+> 0) = —8n(n’ — 5n + 4)ds 
Crv(6; 3, 3, 0, --> 0) = —3(3n + 2)d, 
Crw(6; 3, 2, 1, 0, «+= 0) = —3(n* — 3n? + 4n — 2)ds 
C1w(6; 3, 1, 1, 1,0, «+> 0) = 4(3n® — 11n? + 12n — 4)éd, 
Crw(6; 2, 2, 2, 0, «+= 0) = 2(n* — n® + Gn? — 14n + 8)ds 
Crw(6; 2, 2, 1, 1, 0, ++ 0) = —3(n* — 7n® + 24n? — 38n + 20)d, 
C1w(6; 2, 1, 1, 1, 1, 0, --- 0) = 3(n* — 10n* + 35n? — 50n + 24)d5 
Crv(6, 1,1, 1,1,1,1,0,--- 0) = —ds(n* — 10n*® + 35n” — 50n 


+ 24)és 


A partial check on these coefficients has been made. This was ac- 
complished by taking f(u) as the normal density function and computing 
the coefficients, D;y . Values of D;y and Cy were substituted in (4.13) 
and the values of these integrals were in turn used on the expressions 
for Gy(t) given in (3.11) and (3.13). The series were checked by expanding 
the “Student” density in powers of {’ and integrating that power series 
term by term with respect to ¢ over the range —« to &, say. This 
should give us Gy(é&) and in fact the first four terms are identical with 
the results obtained by the general method of this paper when ¢ is re- 
placed by &. 
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ILLUSTRATIVE EXAMPLES 


We now consider applications of the series expansion of the cumulative 
distribution function of ¢ to two non-normal populations. They have the 
Cauchy and the “squared hyperbolic secant” densities for which f(w) 
takes the forms 


(6.1) (C) e1+u)2 and (S) i sech’ u/2 


respectively. Standard forms may be used since the distribution of ¢ is 
independent of location and scale parameters. 

We first consider the Cauchy density and require expressions for its 
derivatives. A result of J. Bertrand (1864) expresses the general derivative 
of tan” ‘u in compact form and is useful here since f(u) = x‘ a tan ‘u/du. 
Bertrand’s result is in slight error but, with the proper correction, 
(6.2 C) fu) = wy! (—De sin **’y sin (y + Dy, 
where sin y = (1 + wu)” and cos y = u(1 + wu)”. Upon substituting 
u for p/+/N in (4.11), with the help of (4.9), we have D,y expressed in 
terms of the derivatives of f(u). Then, we use the simple transforma- 
tion, uw = cot y and substitute the values, (6.2 C). It follows that 


9 


a/2 
Diw(k; Qo, +++ Qn) = 2 (— sated | cos”** y cosec y 
0 
(6.3 C) ' 
4 sin (a; + y\ dy. 
i=0 


These integrals are evaluated by expanding the trigonometric functions 
of the multiple angles in terms of sin y and cos y. The values required 
are tabled below. 

Derivatives of the squared hyperbolic secant density may be obtained 
simply by repeated differentiation. Even-order derivatives are poly- 
nomials in sech’ u/2; odd-order derivatives are products of such poly- 
nomials and tanh u/2. Integrals, D,;y , vanish unless the sum of the orders 
of the derivatives is even. Thus, when wu replaces p/+/N in these inte- 
grals, tanh u/2 appears with an even exponent and such functions may 
be expressed in terms of sech’ u/2. It follows that all D,;y depend on in- 
tegrals of the form, 


wo 


(648) I(p,q) = [ x” sech™ 2/2 dr = 2 | x(1 + cosh x)“ dz, 
/0 0 


where p is a positive integer or zero and q is a positive integer. 
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I(p, q) is not in general an elementary function; however the following 
recursion relations are useful: 





2(q-1 
(p,q) = 24—® 1p, q — 1) 
” (2q — 1) 
(6.5 S) 2n( 1) 
<p Pp = ‘ ‘ 
gies Se Ra Bh >2 
(@ — 1)(29 — 1) (p q ) Pp ’ 
and 
2(q — 1 
faa = ae 14-1 
(6.6 S) 7 : 
— 2'P i 2?" (1 + eosh x)’ sinh x dx 


The second relation is used when p = 1 and its last term is then easily 


TABLE I 


VALUES OF I(p, q) 








pt 2 3 4 5 6 7 

1 1.18173 | | 

2 0.84963 | | | 

3 1.24309 | 0.41099 | 

4 0. 86625 0.41994 | 

5 | 4.56605 | 1.54594 0.91750 | 

6 4.13962 2.23589 1.57450 

7 29.25208 6. 80901 2.44526 1.82169 

8 17.08400 8. 28641 2.98001 0.48996 
9 11.91614 4.43072 | 2.57293 


integrable. Repeated use of the reduction formulae will terminate with 
one or more of the integrals’: 

< 

| z"(1 + cosh x) dx = 2(2"* — 1)r"B,, 

Jo 


x.) 


| z"**(1 + cosh x) dx 
0 


me ots Be ee 
(6.7 S) = 41 (2a + 2) ara =} be (2; rs 1) (2a of a 


j=0 


IIV 


1, 


6 Bg is the at” Bernoulli number, and, of the various definitions of these numbers, the one used here 
is that Ba is the coefficient of 22*/(2a)! in the expansion of re*/(e7 — 1). Bi = 1/6, Bz = 1/30, Bs = 1/42, 
Bs = 1/30, Bs = 5/66. 
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[ «(1 + cosh x) dx = 2 log, 2, and 
0 


[ (1 + cosh x)~ dx = 2’6(b, b). 
0 


The values of I(p, g), which we shall require, are shown in Table I. 

The coefficients, D,y , are given for both the Cauchy and the squared 
hyperbolic secant densities in order in (6.9 CS). We abbreviate by 
setting 








(68S) 8 (? z, +t3) = b,j) and I(n+i,n+j) = dl, 9. 
Then we list: 
(69CS)  k=0 w= N™e*, Ay = NPQ, 
Diw(0; 0, «- 0) = yw b(1, 1)/2; And, 1). 
k= 2 y= NOMA Ay = NOM 0 
Dix(2; 2, 0, +++ 0) = 2{3b(3, 1) — 2b(3, 3)}; 
As{d(2, 1) — $d(2, 2)}. 
Diw(25 1, 1, 0, + 0) = 2y B(5, 1); 2Aefd(2, 1) — a(2, 2)}. 
k=4 y= NOWARN Ay = NOH Pg p44 
Diw(4; 4, 0, ++ 0) = yf $b(5, 1) — 100(5, 3) + 86(5, 5)}; 
Aufd(4, 1) — 4Pd(4, 2) + 4PdC4, 3)}. 
Din(4; 3, 1,0, +++ 0) = yf86(9, 1) — 40(7, 1)}; Av 4a (4, 1) 
— 16d(4, 2) + 12a(4, 3)}. 
Din(43 2, 2, 0, «+ ,0) = yel$0(5, 1) — 120(5, 3) + 80(5, 5)}; 
Au{6d(4, 1) — 18d(4, 2) + Fd(4, 3)}. 
Din(4; 2, 1, 1,0, +++ ,0) = y«{60(7, 1) — 80(7, 3)}; 


A,{ 12d(4, 1) — 30d(4, 2) + 18d(4,3)}. 


Diw(4; 1, 1, 1, 1, 0, stats , 0) - 8y.b(9, 1); A,{24d(4, 1) 


— 48d(4, 2) + 24d(4, 3)}. 


k = 6 Ye = ay, As = rrr fe 











D 


L 


D 
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Dyyw(6; 6, 0, alleles 


D,w(6; 5, 1, 0, «+ 


Diw(6; 4, 2; 0, sl 


Drx(6; 4, 5, ‘. 0, owe 


Diw(6; 3, 3, 0, wt 


Dyw(6; 3, 2, 1, 0, +++ 


Diw(6; 3, 1, 1, 1, 0, ++ 


D,x(6; 2, 2, 2, 0, ++ 


t AND F STATISTICS 


,0) = y«{32b(13, 1) — 400(11, 1) 
+ 12b(9, 1) — 3b(7, 1)}; 
Ae{d(6, 1) — §3d(6, 2) + 105d(6, 3) 
— 2454(6, 4). 
, 0) = ye{32b(13, 1) — 32b(11, 1) + 6b(9, 1)}; 
Ae{6d(6, 1) — 96d(6, 2) + 225d(6, 3) 
— 135d(6, 4)}. 
0) = yel*Pb(7, 1) — 40b(7, 3) 
+ 64b(7, 5) — 32b(7, 7)}; 
Ae{15d(6, 1) — 135d(6, 2) + 
11254(6, 3) — ®45d(6, 4)}. 
, 0) = ye{10b(9, 1) — 406(9, 3) + 32b(9, 5)}; 
Ae{30d(6, 1) — 255d(6, 2) 
+ 450d(6, 3) — 225d(6, 4)}. 
, 0) = ye{32b(13, 1) — 326(11, 1) 
+ 8b(9, 1)}; 
Ae{20d(6, 1) — 140d(6, 2) 
+ 300d(6, 3) — 180d(6, 4)}. 
0) = ye{24b(11, 1) — 1209, 1) 
— 32b(11, 3) + 160(9, 3)}; 
As{60d(6, 1) — 330d(6, 2) 
+ 540d(6, 3) — 270d(6, 4)}. 
,0) = ye{32b(13, 1) — 16b(11, 1)}; 
As{120d(6, 1) — 600d(6, 2) 
+ 840d(6, 3) — 360d(6, 4)}. 
,0) = ye *,'b(7, 1) — 546(7, 3) 
+ 72b(7, 5) — 326(7, 7)}; 
Agf90d(6, 1) — 405d(6, 2 
+ 349}5d(6, 3) — 2925d(6, 4)}. 


19 











20 THE VIRGINIA JOURNAL OF SCIENCE [January 


Diw(6; 2, 2, 1, 1, 0, «++ ,0) = ye{18b(9, 1) — 48b(9, 3) + 326(9, 5)}; 
As{180d(6, 1) — 720d(6, 2) 
+ 945d(6,3) — 405d(6, 4)}. 
Drw(6; 2,1, 1,1,1,0,--+ ,0) = ve{24b(11, 1) — 32b(11, 3)}; 
As{360d(6, 1) — 1260d(6, 2) 
+ 1440d(6, 3) — 540d(6, 4)}. 
Diw(6; 1, 1, 1, 1, 1, 1,0, «++ ,0) = 32y¢b(13, 1); 
A«{720d(6, 1) — 2160d(6, 2) 
+ 2160d(6, 3) — 720d(6, 4)}. 


The distribution functions of 4 may now be obtained. The sets of 
coefficients (6.9 CS) are combined with the set (5.3) in (4.13). For the 
Cauchy density, 


if a a d B(®) = P(2z) II dz; 


o ap b= 


_ Sra No-Diag £ ln + ') ae 


Qa 2 2 





n+1 
n + 2’ 
_ 3(n + 1)(—2n* — 5n® + 29n + 14) 
i (mn + 2)*(n + 4) 
15(n + 1)(—14n® — 375n* — 146n° 
we = a _+ 3927n* + 3274n + 48). 
(n + 2)?2(n + 4)?(m + 6) 


s = 0,2,4,6, w=1l m= 


(6.10 C) 


W4 


In the second example, for the squared hyperbolic secant density, 


| te | d B(®) P(2) Il ae sl Sn-1 NOP2y, | 
— 20 — oo | d=0 1 


de 2 


s = 0, 2, 4, 6, ho = (n,n + 1), 


—(n + 1) 
2 


Ae = 


I(n + 2,n + 2), 
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3(n + 1) 2 , 


—2nI(n + 4,n + 2)}, 


de = M+ patrin + 6,0 + 2) 


23(n + 2)(n + 4) 
(6n* + 68n” + 24n — 8)I(n + 6, n + 3) 

+ (n‘ + 17n*® + 104n* — 230n + 288)I(n + 6,n + 4)}. 
The cumulative distribution function of ¢ for the Cauchy distribution 
is at once found by substituting the values of (6.10 C) in (3.11). S,_, 
is replaced by its value (5.1). When ¢ > 0, 


n+1n+1 
NFP ng??? B ( 2 ’ 2 ) 


(N+) /2 _[n 
I (3)e 
fy mety G) 4 mn + (14 + 29n — Sn — 2n’) 
(6.12 C) 2(n + 2)?\# 272'(n + 2)2(n +4)? — 
n'(n + 1)(48 + 3274n + 3927n° 
(i) + = 1460n* — 375n* — 14n’) 
2°3 !(n + 2)?(n + 4)?(n + 6)? 


Gy(t) = 1 = 











Similarly, for the squared hyperbolic secant density, 


y edad ial +2)/ =n! n—2)/2 


Gu) = 1 — —sar@ae 
! (n,n+1)_ n'I(n + 2,n + 2) S) 
(n + 1) 2.2!(n + 2) t? 





4 Bnit(n' + Gn + 2)I(n+ 4, n+ 3) — 2nl(n + 4,n + 2)} 


274!(n + 2)(n + 4) 
(€ 13 8) 





; (x) - 15n* . 
i 2°6!(n + 2)(n + 4)(n + 6) 
{4n7I(n + 6, n + 2) — (6n* + 68n" + 24n — 8) 
-I(n + 6,n + 3) + (n' + 17n* + 104n° — 230n + 288) 


-I(n + 6,n + 4)} (x) a a) |. 
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Gx(t) has been evaluated for sample sizes up to five. For the Cauchy 
density, we have the series, 





; 2 { 1 1 
(ft) =1— = 0 
Go(t) = 1 wi\! + ge + 5am t ase + 0(t }, 
le 3/3 f 3 3 29 - 
G;(t) =l1- 7p \! +3 at 16 73 + D560 + OCF , 
(6.14 C) ( mn . a 
16 27 2332 -s 
14() =1-—- a O(E 
Gu) = 1 — Saal! + sop + sega — Bare t ( , 
7 5/5 f 65 355 an 
i(t) =1— —— — ( 
GM) = 1 — soe |! +4 a aae ~ gare tO) 


and for the square hyperbolic secant density, we have 
_ 0.29543 f, _ 0.11983 


G2(t) = ; = 
0.01686 , 0.00233 ? 
+ + + ol }, 
G.() = 1 — 043342 {1 _ 0.76466 
e e 
(6.158) ‘ = » — + oe}, 
0.55910 4.018: 
or -_ {1 me . 20 


2.14852 , 3.39309 -s 
++ + O(t \, 

t f° 
1.81014 { 12.49756 - 


a <-i- =| " ———— + Ott a. 





The values of [(p, q) in Table I have been substituted in (6.13 5) to 
obtain this second set of distribution functions. 

The probability that | ¢| exceeds some & may be obtained directly 
from Gy(t) for both examples since each has a symmetric density func- 
tion.’ Figures I and II show the probability of | ¢| exceeding t plotted 
against %& for the Cauchy examples of (6.14 C). A logarithmic scale is 
used for é& and values for the “Student” distribution are also shown 
for the purpose of comparison. Figures III and IV show similar plots 
for the squared hyperbolic secant example. 


7 For populations with skew density functions it is necessary to use (3.13) as well as (3.11). 
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THE GENERALIZATION TO TWO OR MORE SAMPLES 


We have obtained general expressions for the density function and the 
cumulative distribution function of ¢ for a sample of size N from popu- 
lations whose density functions f(w), satisfy certain: specified con- 
ditions. The methods employed may be generalized to include the 
analogous two-sample and k-sample problems when the populations 
are similarly described. The sequences of transformations which have 
been used permit the required generalizations and, when more than one 
sample is considered, the linear and spherical transformations are per- 
formed in the same order as before.. We state the results obtained and 
again refer the reader to footnote 2. 

In this section, = denotes a summation over 7, with 7 taking the values 
1, --- ,k and S denotes a summation over a for the values 1, --- , N; 
II and M represent respectively the corresponding products and 2’, 
S’, Ill’ and M’ are sum and product symbols defined as above with the 
exception that their upper limits of operation are all reduced by unity. 
N; is the number of observations in the 7‘" sample. 

The ¢ statistic used to test hypotheses on the means of two samples, 
Ni, *** , ty,i, 7 = 1, 2, when the samples are assumed to come from. 
populations whose density functions differ at most in location, is 


(7.1) t = N~?(nN, N2)"*(%, — %)/{2S(tei — ¥)7}"”, 
where N = N, + Nz, n = N — 2, and © is the described summation 
with k = 2. On the hypothesis that the means of the populations from 


which the samples are drawn are equal, each population may be de- 
scribed by the density function, f(u). Then, the density function of ¢ is 


G'(t) = N(NiN,)"”” si ail | ve { 


(7.2) [TIM 'f{u — (—1)*(Ni/Ne)*?P(u) + Vn fUai}] 
-[Mf{u — (—1)'(Ni/N2)*? P(u) — nt S'uas}] 
P(u) du 11M’ dua; 
with 
P(u) = N~"*{3S8'u2; + 3(8’uai)"}™”. 


The cumulative distribution function of ¢ may be expanded in a series 
as before. This expansion is 
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n 


yr \1/2 n/2 
Gi) = 1 — N(NiN2) T(n)n be at Qs)! 
(n+ 28-1) mien f(a BO) , \ 
( Qs )e [ . a lenol (™) GulLM duai 


when ¢ > 0 and provided the series converges. Terms in the series may be 
simplified since 


(7.3) 


[- f a B(®) | P(u) dullM’ duai = 1(NN,N2) (2s)! 


dbs e 
74) Ye (™ + Or — 1 Ne + 28 — Or — ') 2 
= 2 2 @o+-+-+0n,—1=2r 
2» Ciw ,(2r; a)Crw,(28 — 2r; b) Diw,w,(2r, 2s — 27;a,b) 


bot: --+bn,,— 1=28—2r 


with 


foo) @o 
N+2s—2 
Drw,wn,(2r, 2s — 2r; a,b) = | [ F a 
a " 


(7.5) | f*-"(u + ica HY 


' 
a Qa-1 








f (B- 1) a T 1/2 
[ut {u (Ni/NN2) Ma dp. 
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The coefficients C;y, and C;y, are identical with those in (4.12) when 
N is replaced by N, and N,. When ¢ < 0, a similar series is obtained. 
Then, we define Dfy,v, like Dry,v, with the exception that —p replaces 
p in the arguments of f, and the cumulative distribution function of ¢ 
becomes 1 — G*(|¢|) where G*(|¢|) is G(é) defined by (7.3) and 
(7.4) with |¢| replacing ¢ and Div, replacing D;y,~,. The integrals 
D.w,n, and Dfy,v, are of multiplicity two and must be evaluated for 
each density, f(u), separately. 
The analogous k-sample test. depends on the statistic 
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where N = =N;,n, = k —:1, and n, = N — k. The density function of 
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before. The explanation of (7.9) is complete when we write 
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D, is a k-fold integral which must be evaluated for each population 
considered. It depends indirectly on d; and s since fd; = s and i 
Gai = 2d; ; a; represents the vector, (do; , +++ , @n,-1,:). 

It may be shown that when n, = 1, the expression (7.8) for G(F) 
is equivalent to the cumulative distribution function of @ obtained by 
combining G(t) in (7.3) with the similar expression for the cumulative 
distribution function of ¢ when ¢ < 0. 


DISCUSSION 


The examples on the Cauchy and the squared hyperbolic secant 
populations indicate that as the sample size increases, more terms of the 
series may be required to calculate the probabilities that specified values 
of the statistics be exceeded. The rapidity of convergence appears to 
decrease as the sample size increases and the value of the statistic de- 
creases. Due to the complexity of the methods, a consideration of the 
convergence of the distribution functions was unsuccessful. It is hoped 
that further work may be forthcoming on this point. 

It is interesting to note that the coefficients, C;, , appear in all three 
problems and need only to be evaluated once since they are independent 
of the population considered. These coefficients have been evaluated 
for the first four terms of the series but become increasingly numerous 
and complex as more terms are considered. Additional values should be 
obtained but require improved computing methods. 

The series expansions may be used to obtain the power functions of 
the statistics ¢ and F. f(u) may be replaced by f(u — a), where a is a 
measure of location, and the coefficients C;y will be unchanged. When 
both scale and location parameters are present, the measure of non- 
centrality will be a/c, where c is the scale parameter. Some difficulty 
may be introduced in the evaluation of the D-integrals but their multi- 
plicities will be unchanged. 
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SUMMARY 


The probability density functions of ¢ and of F have been expressed 
as multiple integrals of multiplicity one less than the sample size and 
whose ranges of integration are unlimited. The methods are appropriate 
for non-normal density functions satisfying rather mild conditions. 

The Taylor expansion is used to expand the cumulative distribution 
functions of ¢ and of F in power series in ¢' and F‘. The series are 
based on the density functions described immediately above. The coef- 
ficients in the expansion depend on integrals whose multiplicities are 
equal to the number of groups of observations in the sample. Two 
examples are included in section 6. They substantiate the opinion that 
the effects of moderate departures from normality may not be serious in 
the analysis of variance. 

Certain problems associated with this work ‘still require solution. 
However, this is the first study by analytical methods of the effects of 
non-normality for a general class of populations using specified density 
functions. 
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EFFECTS OF INSULIN IN VERTEBRATES 


A Comparison of the Gross Effects of Insulin 
Injection in Normal Mammals, Birds and 
Poikilotherms 


D. R. H. Gourtey 
Department of Pharmacology 
University of Virginia Medical School 


In the period which has elapsed since Banting and Best (1922) 
first demonstrated the internal secretion of a hormone, insulin, by the 
islets of Langerhans in the pancreas, a vast amount of research has been 
directed at elucidating the various functions performed by the hormone 
in the living organism. Much of this research has necessarily involved 
the use of animals or specific tissues of animals of a variety of species. 
The results obtained have been frequently applied a priori to the human 
species. That different species react differently to removal of the pan- 
creas has been pointed out by Peters and Van Slyke (1946). It is the 
purpose of the present communication to review the effects produced by 
administration of insulin to intact animals of different species. 

The concentration of insulin in the pancreas is reasonably con- 
stant from species to species. Normal values have been reported for fish, 
chickens, ducks, mice, rats, guinea pigs, rabbits, cats, dogs, cattle, pigs, 
sheep, horses, monkeys, chimpanzees and men (Haist, 1944). The in- 
sulin preparation used almost exclusively at the present time is isolated 
from beef and pork pancreas. While this hormone preparation has quali- 
tatively similar effects in all vertebrates, it has been shown that there are 
gross quantitative differences in the reactions of various species to in- 
jection of insulin (Macleod, 1923, 1926). 


Effects in mammals.—The main effects of injecting insulin into 
mammals are very striking. The blood glucose level falls rapidly and the 
animal is said to be in a state of hypoglycemia. If the level of blood 
glucose falls below a definite value, which differs from one mammal to 
another, highly characteristic symptoms supervene. These are known as 
the “hypoglycemic symptoms” and, in general, are characterized by con- 
vulsions or coma, or both. 

In the rabbit, the fall in blood glucose sets in almost immediately 
after the injection of insulin so that after thirty to forty-five minutes it 
stands at about one-half its normal value. Definite convulsions appear if 
the blood glucose level falls below this value (45 mg. %) which usually 
occurs following a dose of one to five units’ per kg. body weight. The 


1. By definition 1 mg. of zinc insulin crystals contains 22 international units 
of insulin. 
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typical effect of insulin on the blood glucose of the dog is very similar 
to that observed in rabbits, with the difference being that the fall is 
pet as steep. The minimum blood glucose level (about one-half the 
normal lieve!) is reached in two hours or more, and when the dose is 
about 23 units per kg., convulsions are observed in from four to five 
hours. The cat and rat react in much the same manner as the dog. In 
man the blood glucose decreases more slowly after insulin injection. When 
the concentration of blood drops to 70 or 60 mg. %, the general symp- 
toms may become so severe that it is necessary to administer glucose as 
an antidote. However, in man, variations between individuals are very 
great in this respect. A blood glucose level of 35 mg. % is usually 
accompanied by unconsciousness. Among the mammals on which good 
information is available, the mouse, rat, dog, and man display about 
equal sensitivity to insulin while the rabbit is three to four times as 
susceptible. 


Effects in birds—It has been established that birds, as a whole, are 
affected less by insulin than mammals (Chen et al., 1945, Sordelli et al., 
1923). However, even among different species of birds, there is a wide 
difference in sensitivity. Both the horned owl (Bubo virginianus) 
(Nelson et al., 1942) and the Pekin duck (Chen et al., 1945) are rela- 
tively sensitive to insulin, the latter having convulsions as early as two 
hours after receiving adcquate doses of insulin. The pigeon is more 
tolerant to insulin than the duck, and some specimens can withstand up 
to 150 units within 24 hours, showing only a 50% reduction in blood 
glucose (Stare and Baumann, 1940). Domestic fowl (white Leghorns, 
Plymouth rocks, Rhode Island reds), however, appear to be very resistant 
to insulin; doses of 10 to 130 units per kg. may lower the blood glucose 
level to 25% of the normal, but although the birds get weak when the 
blood glucose is so low, no convulsions appear (Cassidy et al., 1926). The 
observations of Chen and his co-workers (1945) on roosters (white Leg- 
horn) are in accord with this finding. They report that, of nine roosters 
injected with doses varying from 100 to 500 units of insulin per kg., 
only one had ataxia and two fell on their sides, but all appeared normal 
the next day. 


Effects in poikilotherms.—Cold-blooded vertebrates also exhibit a 
marked resistance to insulin. Moreover, it was noted very early that in 
cold-blooded animals it was impossible to correlate the symptoms usually 
associated with hypoglycemia with the blood glucose level. In some 
species the blood glucose level may fall, but convulsions do not appear 
until much later, while in other species insulin injection may cause con- 
vulsions but the fall in the blood glucose, if it occurs at all, is not 
detected for a matter of days later. 


Houssay and his co-workers (1924, Sordelli e¢ al., 1928) found 
that injection of large amounts of crude insulin into a variety of cold- 
blooded animals of the Argentine (including some reptiles) eventually 
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lowered the blood glucose level and produced convulsions. But the 
appearance of the two effects did not coincide in these animais, and ail 
animals died several days after the injections. Small crocodiles (Caiman 
sclerops), for example, were hypoglycemic one day after crude insulin 
injections, showed hypoglycemic symptoms on the third or fourth day, 
and died in seven to ten days (Houssay and Rietti, 1924). 


The cold-blooded animal which has been most thoroughly studied 
in relation to insulin effects has been the frog. Abderhalden (1924) was 
one of the first to note that frogs tolerate remarkably high doses of 
insulin. In more detailed studies, Olmsted (1924, 1927) and Huxley and 
Fulton (1924) observed that frogs at room temperature did not convulse 
earlier than 24 hours after insulin injection. The insulin convulsions in 
frogs have been described accurately, and it has been observed that the 
onset of the symptoms is greatly accelerated by warming the injected 
frogs (Huxley and Fulton, 1924, Olmsted, 1924). The same result ob- 
tains when insulin-injected toads are warmed (Olmsted, 1926). More 
detailed studies of the effects of insulin on the blood glucose level and 
the onset of convulsions were made by Gabbe (1925) with the European 
frog (Rana temporaria) and by Barlow, Vigor and Peck (1931) with 
“grass frogs” (probably Rana pipiens). In both species it was found 
that the appearance of convulsions anteceded, and was apparently inde- 
pendent of, the drop in blood glucose concentration. Gourley and Fisher 
(1951) confirmed these observations in Rana pipiens, and found that the 
injection of 38.0 x 10~ units of insulin per ml. of body fluid caused 
convulsions to occur within 24 hours and a 50% decrease in the blood 
glucose concentration after 48 hours. A concentration of 10~ units of 
insulin per ml. of body fluid had none of these effects. Olmsted (1924), 
who thought that the blood glucose level of insulin-injected frogs re- 
mained normal, failed to observe hypoglycemia only because he made 
his determinations at the time of convulsions. It will be noted that the 
order of appearance of the two effects in mammals is reversed in the frog. 


Olmsted (1924) described peculiar symptoms, considered by him 
to be related to hypoglycemia, which were observed two days after the 
injection of insulin into fresh water catfish (Ameiurus nebulosus). 
McCormick and Macleod (1925) observed similar symptoms in the sculpin 
(Myozocephalus) and the sea raven (Hemitripterus) several days after 
insulin injection but were unable to demonstrate that the insulin injection 
caused the normally low blood glucose to become further lowered. Macleod 
(1926) was inclined to feel that asphyxia rather than hypoglycemia may 
have been partly responsible for the symptoms observed. The possibility 
of asphyxia was eliminated in the carefully controlled investigation of 
Gray and Hall (1930) which indicated that the very active fishes, 
menhaden (Brevoortia tyrannus), common mackerel (Scomber scombrus) 
and bull’s eye mackerel (Pneumatophorus colias), suffered insulin shock 
in a much shorter time than did the moderately active scup (Stenotomus 
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chrysops). The more sluggish bottom feeders, sea robin (Prionotus caro- 
linus), puffer (Spheroides maculatus) and toadfish (Opsanus tau), 
showed no external evidences of the effects of the hormone. Insulin 
appeared to reduce the blood glucose concentration of fishes in much 
the same manner as in mammals, except that a longer time was required 
for the effect to occur. It is curious that in fishes, as in birds, the blood 
glucose may be lowered without any signs of convulsions appearing. This 
fact led Gray and Hall to conclude that insulin shock in fishes does not 
have so much significance as has been attributed to insulin convulsions in 
mammals. 


Summary and conclusions.—The association of the blood glucose 
level to the incidence of hypoglycemic symptoms in mammals suggests a 
causal relationship between the two (Himwich, 1944). The suggestion 
appears to be sound in view of the fact that the hypoglycemia which 
follows removal of the liver is also accompanied by the same symptoms. 
The symptoms are relieved immediately in both instances by the admin- 
istration of glucose. The facts reviewed above, however, indicate that the 
same relationship does not hold in birds and in poikilotherms. Even in 
mammals there are many exceptions (Barlow et al., 1931, p. 236; Gabbe, 
1925, p. 215; Harris, 1944, p. 354; Peters and Van Slyke, 1946, p. 229). 


It is clear then that distinctions between the reactions of different 
species to insulin exist. Mammals appear to be the most sensitive to 
insulin while birds are most resistant to the hormone. An extensive study 
of the effect of insulin in reptiles has not yet been made but an amphibian, 
the frog, appears to be somewhere between mammals and birds in insulin 
sensitivity. Peters and Van Slyke (1946) suggest that these differences 
may be related to differences in the inherent metabolic processes upon 
which the action of insulin is superimposed. These processes are the 
result of evolutionary adaptations and insulin does not bring them into 
existence but may, by influencing the rates of certain ones, alter the 
direction of the metabolism as a whole. 


It is interesting to speculate on the possible reason for the species 
difference in reactions to insulin. The birds of the present period consti- 
tute an extraordinarily successful and highly diversified group, probably 
at the height of their evolutionary career (Newman, 1939). It seems 
conceivable, therefore, that birds might have developed an independence 
of whatever reactions are controlled by insulin. In mammals, on the 
other hand, certain processes might conceivably have become more de- 
pendent on the action of this hormone. A similar evolutionary change has 
been suggested by Beadle (1945) to account for the manner in which 
successive generations of microorganisms lose some of their synthetic 
abilities and become more exacting in their nutritional requirements. The 
ancestry of both birds and mammals has been traced back to primitive 
reptiles, which in turn developed from the primitive amphibians. As has 
already been noted, frogs are more sensitive to insulin than birds but 
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less sensitive than mammals. It is necessary to suppose in the present 
speculation that, in respect to its insulin sensitivity at least, the modern 
frog has not changed greatly from the primitive amphibians from which 
it descended. 

Whatever reasons lie behind the species distinctions in the gross 
effects of insulin, it is clear that in animal experiments designed to 
elucidate the mechanisms of the hormone’s action analogies to other 
species must be applied with caution. 
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The Exoskeleton and Musculature of the 
Appendages of the First Three Abdominal 
Segments of Cambarus Longulus Longulus 
Girard (Decapoda, Astacidae)’ 


C. W. Hart, Jr. 
University of Virginia and Florida State University 
INTRODUCTION 


Since Schmidt (1915) extended the pioneer work of Huxley (1880) 
by rendering a complete account of the musculature of Astacus fluviatilis 
(Astacinae) no comparable study has appeared which deals with this or 
other groups of crayfishes. Since the Cambarinae is a group, restricted 
to the North American continent east of the Rocky Mountains (except 
for introductions by man into California, Japan and Western Europe) 
it seems desirable to examine the musculature of the appendage of this 
related but anatomically distinct subfamily as exhibited by Cambarus 
longulus longulus Girard.” 


The anatomy of the first three abdominal segments, their appendages, 
and a study of the appendicular musculature of this crayfish has been 


found, except for minor details, to be identical with Schmidt’s descrip- 
tion of Astacus fluviatilis. 


The muscular terminology employed by Schmidt is adopted in this 
study. 


I wish to express my thanks to Dr. Horton H. Hobbs, Jr. who has 
read and criticized this manuscript. 


THE EXOSKELETAL CONSTRUCTION OF THE 
ABDOMINAL SEGMENTS 


Several rather distinct regions of the exoskeleton of the abdominal 
segments have been distinguished by previous investigators, and although 
these regions constitute a continuous whole, it is convenient to consider 
each segment to be composed of four basic sclerites—a tergite, a pair of 
pleurites, and a sternite (figures 1-8 and 10). 


The sternum of the first abdominal segment is so constructed as to 
form a heavy transverse ridge which projects ventrally to a marked 
degree. It is on this ridge that the first pair of pleopods articulate 


1. Contribution from the Miller School of Biology, University of Va. 
2. The specimens on which this study is based were collected from Swift Run, 
a tributary of the James River, seven miles west of Stanardsville, Green Co., Va. 


on U.S. Hy. 33, and Keswick Creek, eight miles east of Charlottesville, Albemarle 
Co., Va. on State Hy. 22. 
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(figs. 1 and 10). This sternal region is clearly demarked from the 
pleural sclerites by obvious sutures (fig. 1), and the muscles which con- 
trol the pleopods have their origins on the sternal sclerite. This condition 
is quite different from that found in the other segments of the abdomen 
where the sterno-epimeral juncture is at the mid-point of the articular 
cavities and the muscles have their origins on the tergum (figures 2-6). 
However, the serial homologies of these muscles are quite evident through- 
out the abdominal segments. In spite of the modifications of the first and 
second abdominal appendages, the first three show similarities in their 
structure. 


In view of the fact that the third pleopod is a more generalized 
appendage than the first and second, it is used here as a basis for com- 
parison of the three. 


THE THIRD PLEOPOD 


Gross anatomMy—This appendage, which is attached at the articu- 
lar cavity situated between the sternal and the pleural sclerites of the 
third abdominal segment, consists of a stalk which is composed of a 
short basal joint, the coropodite, followed by a long cylindrical second 
joint, the basipodite (figures 3 and 4). The coxopodite forms a semicircu- 
lar plate, and this together with the basipodite is called the protopodite. 
At the distal end of the basipodite there are borne two narrow, flattened, 
multi-articulate rami, one of which is called the endopodite, and the 
other, which is fixed slightly more distad on the lateral side of the 
basipodite, the eropodite. 


The endopodite, which is larger than the exopodite, is undivided for 
almost half its length. The multi-articular distal portion, the flagellum, 
is narrow and flattened and is divided into a number of small “rings” 
distinguished from one another by slight constrictions of the exoskeleton. 
The exopodite is a similar structure, the undivided proximal portion 
being shorter and not so wide. The edges of both the endopodite and the 
exopodite are fringed mesially, distally, and laterally with long setae. 

Muscutarture (figures 2-4).—Schmidt (1915:238-241) names these 
muscles and gives an account of their function as he found them in 
Astacus fluviatilis. However, since his interpretations of the actions of 
a few of these do not seem to obtain in Cambarus longulus longulus, the 
muscles of this crayfish are described below and their probable functions 
indicated. The names of the muscles used in this paper are the same as 
those used by Schmidt. 


Musculus remotor III. pedis spurii (rem.) has its origin in two 
heads situated somewhat caudad on the broad surface of the tergum. The 
two origins Jie in a longitudina! axis, and the fibers, uniting proximally, 
are inserted caudally on the proximal margin of the coxopodite. It is this 
single remotor muscle which binds the coxopodite to the body and moves 
it backward in the longitudinal axis. 
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Musculus rotator dorsalis basipoditis III. pedis spurii (M. rot.d.) 
originates cephalad to the two heads of the remotor muscle on the broad 
surface of the tergum and is inserted on the proximal preaxial margin 
of the basipodite. It is subequal in size to that of the remotor muscle and 
moves the basipodite laterally. 


Musculus rotator ventralis basipoditis III. pedis spurii (M.rot.v.) 
originates on the sternal sclerite immediately caudad to the articular 
cavity, and is inserted on the cephalic margin of the basipodite mesiad 
to the insertion of M. rotator dorsalis. It moves the basipodite mesially. 
A forward movement is obtained by a simultaneous contraction of the 
two rotator muscles. These two, together with the remotor, can move 
the basipodite in either a transverse or a longitudinal plane, and there- 
fore are able to cause it to rotate. 


Musculus reductor basipoditis III. pedis spurii (M.red.b.) originates 
caudally on the coxopodite distal to the insertion of the remotor, and ex- 
tending through the entire length of the basipodite, is inserted postaxially 
on the distal end of the basipodite. The function of this muscle is de- 
pendent upon the relative positions of the coxopodite and the basipodite. 
This being one of the strongest muscles, it forms, together with the re- 
motor, an antagonist for both of the rotators of the basipodite. 


Musculus adductor endopoditis III. pedis spurii (M.add.end.) orig- 
inates mesially near the proximal end of the basipodite and is inserted 
mesially on the proximal end of the endopodite. 


Musculus abductor endopoditis III. pedis spurii (M.abd.end.) is 
weaker than the above mentioned muscles. A short muscle, it originates 
mesially in the distal third of the basipodite and terminates in a small 
tendon which is inserted on the proximolateral end of the endopodite. 


Musculus abductor eropoditis III. pedis spurii. (M.abd.ex.) orig- 
inates laterally on the proximal end of the basipodite, and is inserted 
laterally in the short basal part of the exopodite. 

These last three muscles, M. adductor endopoditis, M. abductor 
endopoditis, and M. abductor exopoditis behave in such a manner that 
their actions are more easily understood if they are considered together. 
The proximal portions of the endopodite and the exopodite lie together 
so that the movement of one causes the movement of the other in the 
same direction. In this way the three muscles are able to serve as adduc- 
tors and abductors of both the endopodite and the exopodite. 

M. adductor endopoditis and M. abductor endopoditis serve a 
double purpose as adductors for the endopodite and abductors for the 
exopodite, while M. abductor exopoditis serves as an adductor for the 
exopodite and an abductor for the endopodite. M. abductor endopoditis 
is an abductor of the endopodite only in that it moves the endopodite 
laterally past the median axis toward the exopodite. M. adductor 
endopoditis also serves as an abductor of the endopodite in that it assists 
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M. abductor endopoditis. In the same manner M. abductor ezopoditis 
is an abductor for the exopodite only if it moves the exopodite mesially 
past the median axis. The action of these muscles, just as in the case 
of M. reductor basipoditis, is dependent upon the relative positions of the 
joints involved. 


Musculus flagellaris endopoditis III. pedis spurii (M. flag. end.) 
originates laterally in the proximal portion of the basal part of the 
endopodite just caudad to the insertion of M. adductor endopoditis. It 
traverses this segment to the distal end and gives off many sinews, each 
of which is inserted on a “ring” of the flagellum. Thus, by contractions, 
it is able to rotate the end of the endopodite. 


Musculus flagellaris eropoditis III. pedis spurii (M. flag. ex.) 
originates laterally on the short basal segment of the exopodite and is 
inserted by sinews to the “rings” of the flagellum. The action of these 
two flagellar muscles is similar. 


In the male crayfish, the third through the fifth pairs of pleopods 
are practically identical, and in the females the second pair of pleopods 
are unmodifie] and are similar to the third. 


THE SECOND PLEOPOD 


Gross avatomy (figure 9).—The appendages of the first and second 
abdominal segments of the male crayfish are not only larger in size, but 
they are mo«ified in other respects. The second pair are of the same 
general stru.ture as the third and the basal joint (coxopodite) still 
retains its s::mi-circular structure. Here it is composed of two slightly 
overlapping sclerites instead of one continuous sclerite and they are 
designsted a: coxopodite “a” (cx.a.) and coxopodite “b” (cx.b.). The 
next joint, the basipedite, although larger, still retains the elongate cylin- 
drical structi:re found in the third pair of appendages. Extending from 
the protoporlite, just as in the third pair of appendages, are an exopodite 
and an endopcdite, the former remaining similar to that of the third pleo- 
pod and the iatter exhibiting a striking modification. 


The basal half of the endopodite is cylindrical and is similar to 
the basipodite. On the mesial surface of the endopodite, and about three 
fourths its length from the base, is a subtriangular lamella, the apex of 
which is directed dorsally when the pleopod is held against the thorax. 
The distal fourth of the endopodite is a short segmented flagellum re- 
miniscent of the longer multi-articulate portion of the endopodite of the 
third pleopod. 

The ex»podite of the second pleopod is comparable to that of the 
third except that it is of greatcr size. 

Muscutature.—lIn the second pair of pleopods there is little altera- 
tion in the musculature from that pattern found in the third pair, the 
principle diff-rence occurring in the increase in size of the muscles that 
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pull the pleopod preaxially—particularly M. rotator ventralis (figs. 5, 6). 
THE FIRST PLEOPOD 


Gross anatomy.—(fig. 12) The coxopodite, so easily seen in the 
other pleopods, has all but disappeared in the first. The two sclerites, 
ex.a, and cx.b., seen in the coxopodite of the second pleopod, are greatly 
reduced in size and no longer articulate with one another. Ex.a., (fig. 7) 
the sclerite on which the remotor muscle is inserted, lies at the caudolate- 


ral margin of the basipodite and is attached to it by a thick chitinous 
membrane. 


The penultimate podomere, the basipodite, which is an obvious entity 
in the other pleopods, can be distinguished from the distal one only in 
the second form males.’ Here the basipodite is delineated from the 
endopodite by a faintly perceptible line obliquely traversing the ap- 
pendage (fig. 11). In the first form males this line can not be seen 
(fig. 12). The remaining two thirds of the pleopod, which exhibits the 
most bizarre modifications encountered thus far, is homologous to the en- 
dopodite of the other pleopods, as will be shown in a subsequent paper. 


In the first pleopod there is no movable articulation between the 
basipodite and the endopodite, the whole appendage is rigid and the 
above mentiored oblique line serves only as a reminder that at some time 
in the evolution of the crayfish, here probably existed a movable joint. 


Muscvur ature.—(figs. 7, 8) It would seem that the first pair of 
pleopods of the male crayfish would behave like the others; however, 
their movement is largely confined to the longitudinal axis of the body. 
This rigid anpendage is controlled by three strong muscles which cor- 
respond to M. remotor, M. rotator dorsalis, and M. rotator ventralis of 
the other pleopods. The remotor has two heads as does the correspond- 
ing muscle in the second and third pleopods. One head, originating later- 
ally on the sternal sclerite, unites with the mesial head which arises from 
the caudal side of the sternal sclerite just laterad to the mid-point of 
the articular cavity and together they are inserted on the coxopodite. 


This muscle serves to pull the pleopod preaxially between the bases of 
the fifth pair of pereiopods. 


The ventral rotator muscie has its origin on the caudal side of the 
first abdomital sternite and is inserted mesiad of the dorsal rotator mus- 
cle on the cephalic margin of the basipodite. The dorsal rotator muscle 
has its origin laterally cn the sternal sclerite mesiad of the remotor and 
is inserted or the preaxial margin of the basipodite lateral to the inser- 
tion of the ventral rotator muscle. These two muscles are antagonistic 
to the remotor muscle and swing the pleopod postaxially. These three 


3. “Second form males.”—Those individuals which are in the non-breedin 
stage characterized by absence of corneous processes on the distal portion o 
the first pleopod and by blunt hooks on the ischiopodites. Immature males as 
well as alternating adult instars are typically second form. 
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muscles are the only ones found associated with the first pleopod and 
they are homologous to M. remotor, M. rotator dorsalis, and M. rotator 
ventralis four in the other pleopods. 


THE VENTRAL SUPERFICIAL MUSCLES OF THE CAUDAL 
PORTION OF THE CEPHALOTHORAX AND THE 
FIRST THREE ABDOMINAL SEGMENTS 


Figure 2 


Musculus ventralis superficialis thoraco-abdominalis (v.s. th.a.A-B) 
extend from the thorax to the abdomen and fall into two groups, both 
of which originate mesially on the endophragmal skeletal apodemes of 
the last thoracic segment. The first branch, v.s.th.a.A, which is not shown 
in figure 2, is located ventrally to v.s.th.a.B and is inserted medio-ven- 
trally on the interarticular membrane of the joint. The other muscle, 
v.s.th.a.B., is ionger and more powerful and is inserted on the caudomesial 
portion of the sternum of the first abdominal segment. 

Musculi ventrales superficiales abdominis (v.s.a.1-2-3-4) (Milne 
Edwards: La couche superficielle des muscles flechisseurs). The first pair 
of these muscles (v.s.a.1 and 2) originates on the interarticular mem- 
brane between the first and second abdominal sternites. Their origins are 
mediocaudad of the pleura of the first segment and they are inserted 
laterally nesr the caudal margin of the sternum of the second segment. 
The next pair of muscles (v.s.s.8 and 4) have their origins on the inter- 
articular membrane between the second and third abdominal sternites 
just caudolaterad of each side of the sternite of the second segment and 
they are inserted laterally near the caudal margin of the sternite of the 
third segment. 

These muscles are serially homologous and are present between all 
“f the abdeminal sternites. They are attached below the axes of the 
joints and give rise to flexion of the abdomen as a whole by approxi- 
mating the sterna. 


SUMMARY - CONCLUSIONS 


The segmented abdominal exoskeleton of the crayfish consists of 
four sclerites—a tergite, a pair of pleurites, and a sternite—which are 
not always easily distinguished. The ventrolateral plates which project 
below the plane of the sternum may be considered to be tegal folds, 
reminiscent of the branchiostegites which form the lateral portions of 
the cephalothorax. 





The elements which make up the ventral portion of the first abdom- 
inal segment are not easily homologized with those of the other segments 
for not only are the appendages of this segment highly modified, but 
also the sternal area on whicn they are borne is markedly different. It 
is considered here that the heavy ridge which bears the appendages is 
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homologous with the sternite of the other segments. Here the sternite 
has become longer (in the horizontal axis) wii ie ariicular cavities 
not being comparably displaced, thus leaving a considerably broader 
portion of the sternite lateral to each cavity than is found’ in the other 
abdominal segments. 


The muscles of the segment, two of which in the other segments have 
their origins on the tergum, here originate entirely on the sternite. This 
anomaly may be due in part to a gradual thickening and deepening of 
the sternal sclerite resulting in a continual contact and subsequent fusion 
between the muscles involved and the upper part of the sternite. 


Some of the interpretations of muscle action given by Schmidt appear 
to be incomplete when considered in relation to other muscles controlling 
the same joints and the probable functions of these muscles have been 
pointed out in appropriate places above. 


The insertion of the remotor muscle on the first pleopod has been 
referred to as being on the “protopodite.” The author, however, con- 
siders that as in the other pleopods, this muscle is inserted on a sclerite, 
which though greatly modified, is homologous with coxopodite “a” of the 
second pleopod. A more detailed study of this will be brought forth in 
# subsequent paper. 




















EXOSKELETON AND MUSCULATURE OF CRAYFISH 





Explanation of Figures——1. Ventral view of first and second ab- 
dominal segments. 2. Dorsal view of the first three abdominal seg- 
ments and the last thoracic segment. Dorsal portion of the carapace 
and the terga are removed. Viscera and extensor muscles also removed. 
8 and 4. Mcesiocephalic views of the third abdominal segment of male. 
5 and 6. Mesiocephalic views of the second abdominal segment of the 
male. Distal portion of the endopedite and exopodite not shown. 7 and 
8. Cephalic views of sternum, pleuron, and proximal portion of pleo- 
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pods of the first abdominal segment of male. 9. Cephalic view of the 
second pleopod of male. 10. Caudal view of the sternum and the basal 
portions of the pleopods of the first abdominal segment of adult first 
form male. 11. Mesial view of left first pleopod of male (second 
form). 12. Mesial view of left first pleopod of male (first form). 


a.c. articular cavity 

bp. basipodite 

cp. carapace 

cx. coxopodite 

cX.a. coxopodite “a” 

cx.b. coxopodite “‘b” 

end.- endopodite 

enp.s.a. endophragmal skeletal apodeme 
ex. exopodite 

M. abd. end. Musculus abductor endopoditis 
M. abd. ex. Musculus adductor endopoditis 
M. add. end. Musculus abductor endopoditis 
M. flag. end. Musculus flagellaris endopoditis 
M. flag. ex. Musculus flagellaris exopoditis 
M. red. b. Musculus reductor basipoditis 
M. rot. d. Musculus rotator dorsalis 

M. rot. v. Musculus rotator ventralis 

pl. pleuron 

pleo. pleopod 

rem. Musculus remotor 

st. sternum 

ter. tergum 

v.s.al. & 3. Musculus ventrales superficiales abdominis 
v.s.th.a.B. Musculus ventralis superficialis thoraco-abdominalis 
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ON THE PROPERTIES OF THE MULTIPLE 
COMPARISONS TEST*® 


D. B. DuNncAN 


University of Sydney and Virginia Agricultural Experiment Station 
of the Virginia Polytechnic Institute 


The multiple comparisons test is a procedure which has been recently 
proposed by the author (1951) for testing the differences among several 
treatments in an analysis of variance. The purpose of this paper is to 
give a more analytical discussion of the development and properties of 
this test. 

STATEMENT OF THE PROBLEM 

The problem may be stated in general terms as follows. Given a sample 
x = (%,%2, +++, %,) where 2; is drawn at random from a normal popu- 
lation n(u;, 0), 7 = 1, 2, --- , n, and given an estimate s° of the com- 
mon variance o”, which is independent of x and is distributed as o X05 

2 
where n2, the “degrees of freedom” for s’ is also known; a test is re- 
quired for making a decision to accept one of the hypotheses in a system 
S defined in the following way. Let wa, , ua, , *** , Ma, be any one of the 
n! permutations of the population means. Let 





Ma, Maz » Maz; Ma, » Mas Mag Ma; Mag Mag Majo May, *** » Ma, (1) 


— | 











represent any plan of linking adjacent means together in groups of one, 
two, -:- , m means so that no group is completely enclosed in another, 
varying from a plan with n single-mean groups to one consisting of 
only one group of all n means. Let H be the hypothesis stating that the 
left hand member of each pair of means not bracketed together in such 
a plan is less than the right hand member and that the relative rank- 
ings of each pair of means bracketed together is otherwise unrestricted. 
The system S consists of all distinct hypotheses which can be obtained 
in this way. 

In the simple case n = 2, S contains the three hypotheses, uw, ue ; 
Me, wi; ANd yy, we; that is, uw. < we, we < mw, and yw; is unranked rela- 





1 This and the previous paper (1951) include and extend a thesis (1947) approved for the degree of Ph.D 
at the Iowa State College. 

2 For this and a companion paper [Va. Jour. Sci., 2 (4)], presented at the Annual Meeting of the 
Academy, May 11, 1951, Mr. Duncan received the J. Shelton Horsley Award. 
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tive to uw. In the case n = 3, S contains the 19 hypotheses in table 1. 
The number of hypotheses in the system increases rapidly with n. 

In the simple case n = 2, an ideal answer to the problem is provided 
by the symmetrical or two-sided ¢ test. The test we seek for the case 
n > 3 may be regarded as a generalization in a sense of this well-known 





test. 
TABLE 1 
System of Hypotheses. Case n = 3 

M1 » M2 y Ms Hi Me Ms Hi, M2 Ms Hi Me , Ms 

Lif 3 | ae | a 
Mi , M3, Be M1 Bs Be M1, Bs Be Me M1» Ba 

a | L J 

M2, Mi, Bs M2 Mi M3 Be , Bs Ma M3 1 , Be 

L J _—— | a | 
M2, Ba,» Mi Me M3 M1 
Ms, Mi, Me M3 M1 Me 

a 
Ma, M2, Mi Ms Be fi M1, M2» Ms 

Se ee _ a | 


GEOMETRICAL NOTATION 

In considering this problem it will be useful to represent the various: 
comparisons within a sample x = (2, %2,-°-°: , 22) by a sample point 
z= (% , 2, °*** , Zn1+) in a sample space W,_, with axes defined by 
z: = yi/s, where y; is any set of standardized orthogonal comparisons, 
say yi = (Dejan tj — itins)/ViG + 1), 7 = 1, 2, --- , n — 1. Corre- 
spondingly any set of values for the comparisons among the population 
means pi, M2, *** , wa May be represented by a parameter point 
6 = (6:, 0,---, O,1) in a parameter space Q2,_; with axes defined 
by 6, = :/o, where 


n= (Dus - ims) / VieFD, t= 1,2,---,n—1. 


From these choices of axes it follows that the probability element at 
any sample point z in W,,_, is closely related to that of an F distribution 
with n — 1 and ne degrees of freedom and is given by 


f(z; 0) dz = gin + 0) (1 + _ + erie dz (2) 


T'(3ne) (nor) im Ne 


2 
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ni 
° ° ° 2 2 
where F is the variance ratio, >> (y; — 6) /ms, andn = n — 1. 
i=1 


In the important special case n. — ~, (2) reduces to an n — 1 variate 
normal distribution with means 6, unit variances and zero correlations. 

The subset of Q,_, defined by an element of the system S, say by 
Ha; Haz» Ha, *** Hay will be denoted by w(a a »Az-**, an). The construc- 


——_ — + —{—____ 





(@ Gp, ds +++, On> 
—+—_——4 | 


tion of any test 7 will consist of defining a region R 








T’) corresponding to the subset w(a; a2, a3; --- , a), for all subsets in 
if 


the system. The test 7’ may then be described as the rule: if 








zeR (a az, a3--- , Qn; T) 
4 | 
make the decision D (a, a2, a; --- , @,), which consists of accepting 
ae 4 | 


+ 


the hypothesis: 


H (a, G2, Ag ++ , On):0 € w(ay Gz, 3 *** , Gn). 
=~ | —-+—— | 


In addition to the subsets in S we will also need to consider a class 
of “null” subsets: 


Wo(@i , @z,°** ,Qp)iMa, = May = *** = be, 2 pSn, 


and “critical”? subsets 
w (ay pes OFFS Ap) ?Q, a wo (ay Oe ©" 5 py), 


and their corresponding hypotheses, regions and decisions which will 
be denoted by Hi(a; , a2, --- , a»), Ri(a,, a2, --- , ap; T), and 


Dia, @2,°°-, a), t = 0, E. 


Example 1: Figure 1 shows the parameter subsets and sample space 
regions for the familiar a-level two-sided ¢ test in the case n = 2. 

It will be noted in passing that a distinction is being made between the 
hypothesis H(1, 2) which is accepted when the decision D(1, 2) is made 


and the null hypothesis Ho(1, 2): 4: = we which is the central hypothesis 
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tested by this procedure. It will be agreed that acceptance of H(1, 2) 


which merely states ignorance of the relative rankings of uw; and pe, is 
the more accurate interpretation of the decision made by an agronomist 
say, when he concludes that the o»served means of two varieties are 
“not significantly different”. This is not to say however that the null 
hypothesis H)(1, 2) is not the correct one to test. We are merely taking 
the usual view that insufficient evidence to reject this hypothesis is not 
sufficient evidence to accept it. The distinction is of little practical con- 
sequence in this case n = 2, but has to be made in all other cases n > 2 
to avoid contradictions. 





- - ~~~ ~~ all, 2) ——- 
w(1, 2) : w(2, 1) 
wo(1, 2) 
| 
0 0 


Parameter space Q 


R(1, 2; t) Ri, 2; 0) RZ. 1:4) 


Sample Space W, 


FicurE 1. Parameter Subsets and Sample Space Regions, Case n = 2 


The test boundaries + ¢, in figure | are the usual @-level significant 
values of ¢, that is, they satisfy the relation P[—ta < ta, < ta | Ho(1,2)| 
= a, where ¢,, is distributed as ¢ with nz degrees of freedom. 


Example 2. In the previous paper, Duncan (1951), the invalid pro- 
cedure of applying a two-sided a-level ¢ test jointly to the null hypothesis 
H(i, j) for all nC, combinations of 7 and j was termed an a-level multiple 
t test. Figures 2 and 3 show the parameter subsets and the samole sdace 
regions for the case n = 3 of this test. 

Only six of the parameter subsets are shown explicitly ia figure 2 
and these are bounded by the lines uw; = yu; , (7,7) = (1, 2), (1, 3), (2, 3). 
The remaining 13 subsets are obtained as combinations of the ones 
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illustrated. For examples, w(1, 2 3) = w(1, 2, 3) + w( 2,1, 3); 


eas = w(1, 2,3) + w(2, ; +> w(1, 3, 2); and w(1 a, ?, 3) ™ Qe , 


that is, the sum of all the subsets shown. 
The boundaries of the 19 sample space regions in figure 3 are given by 
the six lines (x; — 2,;)/+/2s = + 1, (i,j) = (1, 2), (1, 3), (2, 3). 


WEIGHTING OF ERRORS AND DETERMINATION 
OF AVERAGE RISK 


Two-sipep ¢ Test or H)(1, 2).—In this section we will restrict 
ourselves mainly to the case n = 3, and will first consider the deter- 
mination of an average weighted risk* for an a-level two-sided ¢ test of 
the null hypothesis Ho(1, 2). The parameter subsets in Q, and sample 
space regions in W, for this test are the same as those shown in figure 1, 
except that the boundaries at the points @ = 0 and z = +4, in that figure 
will now be the lines 6, = 0 and z; = +f, respectively in the two dimen- 
sional spaces. 

For the weighting of errors we will use the weight function w.(@, D) 
as defined in table 2, where a(@,) and 6(@,) are monotone decreasing and 
monotone increasing functions of 6, over the ranges 6, < 0 and 6, > 0 
respectively. The importance of the less conservative type-1-like error 
of deciding that a difference is positive when it is not, is measured by 
a(@,) and the importance of the more conservative type-2-like error of 
failing to decide that a difference is positive when it is, is measured by 
b(6,). Since the errors involved here do not coincide exactly with the 
Neyman and Pearson type 1 and type 2 errors we will refer to them as 
type a and type b errors. 

Any one decision in the two-sided test may be said to commit none, 
either one or both types of errors. For the parameter values 6, < 0 for 
example, if the decision D(1, 2) is made no error is committed, hence the 
weight zero. If the decision D(1, 2) is made a type b mistake is com- 


mitted relative to — 6, , hence the weight b(—@,). If D(2, 1) is made it 
may be considered that two errors are committed. One is the type a 
error of deciding that 6; is positive when it is not, the other is the type b 
error of failing to decide that — 6, is positive as it is; hence the weight 


((4;) )+ b(- 6). 


‘See Wald (1939, 1942, and other publications). 
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Because of symmetry considerations it is appropriate to use the same 
weight functions with respect to —6, as with respect to @,. 





e, 
Ny 7, 
69(3,1,2) ty (3,2,0) 
GW (1,3,2) (2,3 ,1) ®, 
P= 3 
G(1,2,3) Gd (2,1,3) 
“Paks 





Ficure 2. Subsets in 22. Case n = 3 


Using the above weighting function, the risk r2(6@) for this two-sided 
test may be written as the sum of a type a and a type b risk as r,(@) 
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Figure 3. Regions for a Multiple ¢t Test. Case n = 3 


TABLE 2 
Weight function w: (6, D) for two-sided test 





Decisions 
Parameter Values . - -_——- ee : a —_—_ 
D(A, 2) | D(1, 2) D(2, 1) 
a< 0, i.e. wi < pe 0 b(—@) a(6;) + b(—6;) 
@, = 0, i.e. wi = we | a(0) 0 a(0) 


6, > 0, i.e. we < wi | a(—6:) + 6(:) b(6:) | 0 
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+ 1,(0), where 


ra(8) = a(,) f (z; 0) dz, 
R: 21> ta 
r(0) = b(—4) f (z; 0) dz, 6, < 0; 
R321; >—ta 
ra(@) = 2a(0) f(z; 6) dz, r(0) = 0, 6=0: (3) 
R:2z,> ta 
ra(0) = a(—6) f(z; 0) dz, 
Riz, <—ta 
rp(@) = b(6;) f(z; 6) dz, A; > 0. 
R: 2, < ta 


If we now let p.(@) represent an appropriate parameter anticipated 
frequency function which will have the symmetry property po(—é@) 
= p,(@), the average weighted risk A» can also be defined as the sum of 
two parts Ag, and A» , which may be termed a type a and a type b aver- 
age risk. For these we may write 


a 


Ay = | ra(O) po(@) dd; Ax» = i ry(0)po(0) dé. (4) 
Qe Qe 


In terms of the above analysis the motivation for choosing an a-level 
two-sided ¢ test for a given requirement may be based on the following 
two points: 

(1) For all a-level tests, and for all symmetrical weight and parameter 
anticipated frequency functions of the above kind, the two-sided ¢ test 
gives the minimum type a average risk A», , the minimum type b average 
risk Ay and hence the minimum combined average risk A2. 

(2) For all two-sided ¢ tests, and for the particular forms of w.(6@, D) 
and p2(@) involved, the test at the chosen level a minimizes A, with 
respect to changes in a. 

With regard to the first point it will be noted that the minimization 
of Aw is equivalent to the familiar property of uniform maximum power 
relative to all unbiased tests. Observance of the second point is restricted 
in practice by limitations on the range of significance levels for tabled 
values of ¢ and also by insu cient consideration of a(@) and b(@) and 
inadequate knowledge of p2(@). Nonetheless it will be agreed that in 
choosing a significance level the aim should be to pick a value such that 
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if a lower, more conservative one were taken, the increase in A» would 
exceed the decrease in Ao, ; while on the other hand, if a higher, less 
conservative level were taken, the increase in A2, would exceed the 
decrease in Aw. 

MUuttTIPLE ¢ Trest.—We will now consider the determination of an 
average weighted risk for the multiple ¢ test. All of the points arising in 
this are applicable also in the multiple comparisons test with which 
we are ultimately concerned. The weight function to be proposed can be 
defined completely, except for one important modification to be treated 
later, as the reduced sum of three two-sided test average risks in the 
following way: 

First, each of the 19 decisions may be treated as the sum of three sub- 
decisions with respect to each of the three differences between the means 
considered a pair at a time. On weighting the errors accordingly we 
may write 


ws(@, Diz) = 2D wildy, Dds), ij = 12, 13,28.) 
¢ ij 


where \ij = (uj; — wi)/W2e and d;; = (x; — 2x,)/+/2s and D(d;;) con- 
sists of the three decisions:—D(i, j): decide nu; < yu; , D(i, 7): decide pi 


cannot be ranked relative to yu; , D(j, 7): decide nu; < u;. The reduction 
factor 2/3 is introduced in allowance of the fact that these differences 
are correlated in such a way that three may be considered as equivalent 
to only two independent differences. Thence, we will obtain a risk func- 
tion r3(@) which is the reduced sum of three two-sided test risk functions, 


r3(0) = rsa(@) + rx(0) 


2 9 (6) 
“3 >. rea(dij) + ; > ra(d:;), tj = 12, 18, 23. 
ij oO ij 


Next, in considering a parameter anticipated frequency function 
p3(@) it will be appropriate from symmetry considerations to choose one 
which gives the same symmetrical marginal distributions p2(A;;) for 
each of the differences \;;. That is, to have po(Ajz) = po(Ais) = Poles) 
= po(—Aw) = po(—Ais) = po(—As3). Then on integrating the risk over 
the parameter space we obtain an average risk A; which is the sum of 
the average risks of two two-sided tests, 


A; = Aza + A » = 2Ao + 2Ay . (7) 


But for the required modification, we could say that the multiple ¢ 
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test had identical properties with those of two independent two-sided 


i tests and thus had a minimum average risk. ( 
The necessity for the modification is most evident when we consider 


the properties of this procedure as a test of the three-mean null hy- 
pothesis Ho(1, 2, 3), that is, Hou; = po = ws. It is only in the region 
R(1, 2, 3; m.t.) that this hypothesis is not rejected. It can be readily 


shown that 


Plz e R(1, 2, 3; m.t.) | Ho(1, 2, 3)] = 


t 


* (dans) da (8) 


Rik 
where f(q3,n,) is the probability density of the studentized range of 
a sample of size 3 with n2. degrees of freedom. This probability is far 
too small and the complementary probability, a; , of wrongly rejecting 
H(i, 2, 3) is far too large to be consistent with the choice of a as the 
allowable probability, a. , for wrongly rejecting any one of the two 
mean null hypotheses Ho(1, 2), Ho(1, 3) and H)(2, 3). For example, 
in the case nz > © and a = .05, we find that a3; = .1223 which is an 
unduly large figure relative to .05. 

It can be argued* that since Ho(1, 2, 3) is the null hypothesis that 
two independent differences are each zero that a; should be made equal 
to 1 — (1 — a)’ to be consistent with the choice of a for ae. 

A relevant graphical picture can be obtained by constructing the 
critical boundaries in W, for two a-level two-sided ¢ tests of the inde- 
pendent hypotheses yp; = we and pw; + we = 2p; in the case np > © and 
a = .05. These are shown in figure 4. The argument is that the null 
region R(1, 2, 3; m.t.) should be as large as the null region R(1, 2, 3; ¢’) 


in this figure. In very loose terms we can say that the attempt by the 
multiple ¢ test to test three differences, where only two independent 
ones are present, is causing overcrowding at the intersections of the 
respective null regions R(1, 2; t), R(1, 3; t) and R(2, 3; #). 

hcl tio Lod 


Weighting for Overcrowding: The problem of including the above 
considerations in the overall average risk function is quite an elusive 
und complex one and the following approach may or may not need one 
or two minor adjustments in the light of further thought on the subject. 
It does however appear to be essentially correct with respect to the 
main issues involved and to lead to conclusions which are not likely to 
be changed by any of the possible modifications. 


‘See companion paper Duncan (1951) section entitled ‘Consistent significance levels’. 
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First, on the circumference of any circle r? = 2} + z3 in W2 we will 
define an overcrowding factor k(r) as the ratio g(r)/h(r) where g(r) 
is the length of arc contained in the intersection of R(1, 2, 3; ?) and 

a | 


| 


~ 














ty 
feet 


RCI, 2,3; 











Ficure 4. Sample Space Regions for Two Independent a-Level Two-Sided t 
Tests. 


R,(1, 2, 3; m.t.), A(r) is the length of are contained in R, (1, 2, 3; m.t.) 
and R, (1, 2, 3; mt.), as indicated before, is the critical region 
W. — R(1, 2,3; m.t.). Thus we obtain 


k(r) = 8r (cos ) / 12r{ cos = = f Geers F ti 


k(r) = 8r (sin! ie *) / 2m, ile <r< V2, 


and k(r) is zero for all other values of r. 
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Next, for any given sample point z, let ZO be the line joining z to 
the origin O, and AB the line passing through O at right angles to ZO. 
Let A’B’ be a line in Q corresponding to AB in W,. That is, if z2 = Kz 
is the equation for AB then 6. = K®, is the equation for A’B’. Let @B’ 
be the perpendicular from any given parameter point 6 to the line A’B’ 
and put + for the length of 0B’, using + or — depending on whether 
6 is on the same side of A’B’ as z is of AB or the reverse side. 

If the decision at z is D,(1, 2, 3) or Do(1, 2, 3), (that is, if it rejects 
or fails to reject Ho(1, 2, 3)), it may be said to imply that the sign of \ 
is positive or unknown respectively. With this in mind we may adopt a 
weight function for overcrowding as defined in table 3. 


TABLE 3 
Weight function for overcrowding. 





Decision 
Values of X P nipsiilideninaiiiatinay 
Do‘1, 2, 3) Dy (1, 2, 3) 
1A<0 0 k(r)-aQ\) 
A>0 


—k(r)bQ) 0 


In doing this we are taking the attitude for example that if \ is not 
positive and D,(1, 2, 3) is made, an additional type a error is committed, 
the importance of which depends on the degree of overcrowding at z as 
measured by k(r) and on the magnitude of \ as measured by a(A). On 
the other hand, if \ is positive and D,(1, 2, 3) is made, this will be a 
correct decision and will offset some of the type 6 errors in the rest of 
the test. The importance of this is also being taken to depend on the 
overcrowding and on the magnitude of \ as measured by b(A). 
Modification in Average Risk due to Overcrowding.—The calculation 
of the modification in average risk due to the overcrowding discussed 
above can be simplified considerably by assuming a spherically sym- 
metrical parameter anticipated frequency function p;(@). Since p;(@) 
is already restricted by the nature of the problem to be symmetrical 
with respect to each of the three differences involved, and since these 
modifications are only small relative to all of the risks concerned, it 
will be seen that this is not likely to lead to appreciable errors. 
Under this assumption it follows that the required modification in 
average risk, which we will denote by Aj, can be expressed in terms of 
the probability densities on the line z; = 0 and within the interval 
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ta < 2 < V/2t.. It may be written a3 the sum o° a type a component 
Aj, and a type b component Ay deined by 


e 


ya | a(@2) | g(r) f (2; 0) dzp3;(0) dé 
62<0 R 


(10) 


, 
Aw 


II 


na b(02) | g(r) f (z; 0) dzp3(@) dé 
62>0 R 

and R is the given region: z, = 0, ta < 2 < +/2t.. Thence the final 

average risk for the multiple ¢ test is 


A; + As = (2A20 + Ase) + (242 — Ass). (11) 


The main problem at hand is to modify this test so as to reduce the 
type a risk at least back to 2A, with a smaller increase taking place 
in type 6 risk bringing it back towards 2A» . 

Tue Muoutrete Comparisons Txsr.—In the proposed multiple 
comparisons test, the hexagonal boundaries of the null region R(1, 2, 3; 


m.t.) are expanded outwards to the likelihood ratio circular boundary 


> 


r = 2Fon,,a,, Where F2,n,,a, is the a3-level significant value of F for 2 
and nz degrees of freedom, and a; = 1 — (1 — a)’. The regions of the 
new test obtained in this way are illustrated in figure 5. 

The six small regions lying outside the new null region R(1, 2, 3; 


j 


m.c.) but inside the previous one R(1, 2, 3; m.t.) require special mention. 


In each of these we do not make any one of the eighteen non-null deci- 
sions belonging to the system S. However, for any sample point z in 
one of these regions we make the new decision that \ is positive, where 
\ is defined’ as above, relative to z. A decision of this type may be con- 
veniently left in the form that two of the means considered as a group 
are higher or smaller than the third. If z falls in the region at the top of 
R(1, 2, 3; m.c.) for example, the decision would be that mw; and pe as 


a group are higher than w;. A convenient notation for this region is 
R(3, 1, 2). The other five regions of this kind are denoted in a similar 
way. The remaining eighteen non-null regions of this test are identical 
with those of the multiple ¢ test except for the new circular boundaries 
that six of them have with the null region R(1, 2, 3; m.c.). 


5 It may be noted that \ = E(c) where c is as defined in the preceding paper (1951), page 185. 














62 THE VIRGINIA JOURNAL OF SCIENCE (January 


Using a similar approach to determine the average risk for this test, 
it is necessary to introduce a compensation-for-overcrowding factor 
k*(r) defined as the ratio g*(r) h*(r), where 7) is the length of are 
on the circumference of the circle r?> = zi + 23 contained in the intersec- 
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FicureE 5. Regions for Multiple Comparisons Test. Case n = 3 


N.B. The six regions R(3,1,2;m.c.), etc. have been exaggerated for clarity by 
decreasing the size of R(1,2, 3; ;m.c.) 


tion of Ril, 2, 3; m.c.) and R,(1, 2, 3; ¢’) and h*(r) is the length of are 


contained in R,(1, 2, 3; ¢’). 
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For overcrowding, the factor is 


k(r) = 8r (sin _ r) / on, V 2F 2,no,c3 < r < V/ be (12) 


and zero for other values of r. For compensation for overcrowding, the 
factor is 


8r (sin fe 4 
k*(r) = ; . 28 la <7 < V2F2,n2,0; (13) 
8r (sin-* _ *) 


4 


and zero for all other values of r. 

On incorporating both the overcrowding factor k(r) and the com- 
pensation factor k*(r) into the average weighted risk in the same way as 
k(r) was used in the previous test, we conclude that the complete value 
for this test is 





A; = 2Ae + A; — A} = Asa + Ay = (2Aoq + Ase said Aja) 


. (4) 
+ (2A» + Ag, — Ago) 
where 
, I a(0.)P’ (6) p3(0) do, 
2<0 
, ee i b(6) P’(6) do, (15) 
62>0 
PO) == | geile; 0 dz, 
4 Jere 
R’ is the region: 2; = 0, V2F2,n.,0; < 2 < VW2ta, 
AL = 4] al@)P*@p) ao, 
62<0 
De es Sa [ b(6:)P*(6)p3(6) 48, (16) 
6o>0 


Pe) = 5 | o*ile; ® ae 


and R* is the region: z; = 0, ta < 22 < ) | 

The complete evaluation of the average weighted risk depends, of 
course, on further knowledge of the weight functions a(@.) and b(62) and 
-on the parameter function p;(@). However, much can be done without 
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this by expressing 2A2, and 2A», in forms similar to (15) and (16) and 
considering the values of the known parts of the complete integrand 
separately for each value of @. Thus we may write 


2Ax = 4 I hs, a(62)P(6)ps(6) de 
and 

ee | _, b() [L — P @lps(@) do, (17) 
where 

PO = | sle:0) de 


and R is the region: z. > t,. Thence we have 


4 a(6.)[P(6) + P’(@) — P*(6)] ps(@) de 


62<0 


Aa 
and 


An =4 J  b(6.)[1 — P@) — PO) + P*@)] ps(6) do (18) 
2>0 

The fairly close agreement between the average weighted risk for 
this test and that of two independent two-sided ¢ tests can be seen from 
a consideration of the smallness of the differences P’(@) — P*(@) relative 
to P(@) for @ < 0 or to 1 — P(@) for 6 > 0. These are shown for exam- 
ple in table 4 for a few selected values of @ in the special case n = 3, 
a= 05,n—- ~. 

From these figures and a consideration of similar details in other 
cases, it will be seen that the type a average risk A3_ is even smaller 
than 2A,, and at the same time A, is also smaller than 2A.5,. 

On further consideration of the mechanics of the above weighting 
system it appears that by choosing the likelihood ratio boundary for the 
region R(1, 2, 3; m.c.) we have succeeded in maximizing the differences 


P’(@) — P*(@) for all values of 6 and thus in minimizing the average risk 
with respect to tests having the same size a; for the null region R(1, 2, 3). 


Furthermore, since Ay = Ao + A» is the minimum risk for a single 
two-sided ¢ test and since 240, — Aza and 2A4 — Ax are both small 
and appear to be at least roughly equivalent reductions, A; = Aga + Aw» 


should be fairly close to the absolute minimum possible for the re- 
quired test. 








b- 
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TABLE 4 


P values.* Multiple Comparisons Test 
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Case n = 3,a = .05,;m — « 
~~ " a= = 
1 | 
\ 0 +1 +5 
62 } 
1 — P(@) .6772 .6772 .6772 
1.5 — P’ (@) — .0467 — .0283 .0000 
P* (6) -0317 .0195 - 0000 
ro = — .0150 — .0088 .0000 
0 P (6) .0250 .0250 .0250 
P’ (6) 0044 .0027 .0000 
— P* (@) — .0044 — .0027 - 0000 
P (6) .0069 .0069 .0069 
—.5 P’ (@) .0012 - 0007 . 0000 
— P* (@) — .0014 — .0008 -0000 
P’ — P* — .0002 — .0001 . 0000 
6 Evaluated using tables (1942) and (1945). 
TABLE 5 
P values. Multiple Range Test 
Case n = 3,a = 05, n2— © 
A 
\ 0 +1 +5 
62 i 
— P’ (@) — .0478 — .0290 - 0000 
1.5 P* (0) . 0367 .0223 .0000 
Pt = | — .Oll11 — .0067 -0000 
P’ (@) .0047 .0029 . 0000 
0 — P* (6) .0047 .0029 . 0000 
P’ (6) .0013 .0008 . 0000 
—.5 — P* (@) .0014 . 0009 .0000 





Pr’ — P* — .0001 -0001 - 0000 
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For illustrative purposes, a set of figures comparable to those in 
table 4 have been determined for the multiple range test obtained by 
expanding R(1, 2, 3; m.t.) to the hexagonal boundaries’ (x’ — 2”’)/s = 


3.0.00 = 2.91849, where x’ and x” denote the highest and lowest of 
the three variates 2 , 2 , 23, and q3,,,,.0975 is the .0975 level significant 
value of the range’ of three variates in the case n. > ~. These are given 
in table 5. The superiority of the multiple comparisons test follows from 
comparing the differences P’(@) — P*(@) for the two tests. 


MULTIPLE COMPARISONS TEST FOR ALL CASES 


The multiple comparisons test for any value of n may be defined by 
a recursion method for constructing its sample regions in W,_,. 

To construct the regions for an n variate test: (a) Consider the n 
possible combinations of the variates taken n — 1 at a time; (b) Divide 
W,1 into n intersecting systems of regions, each representing an n — | 
variate test for each of the n combinations. The regions formed by the 
intersections of these n systems are the required regions except for one 
modification which follows; (c) Adjust the null region R(1, 2, --- , n) to 

oud 








have the boundary: =z; = (n — 1)Fn-1,n.,0, Where the latter’ is the 
a,-level significant value of F with n — 1 and ne degrees of freedom 


“cc 


respectively and a, is the “a-per-degree-of-freedom” level 1 — (1—«)"™". 


The naming of, and thus the decision which goes with, each region 
formed in this way, will follow readily from the regions for which it is 
the intersection, except for certain new regions such as R(3, 1, 2) in the 
‘ase n = 3. These are the regions which fall inside the unadjusted null 
region R(1, 2,---, n) but outside its adjusted and final form 








R(1, 2, --- ,n; m.c.). When a sample point z falls in one of these regions, 








we decide to accept the hypothesis H: \ > 0 where X is defined relative 
to z as in the section on weighting for overcrowding. 

The multiple comparisons test is identical with the two-sided ¢ test 
in the case n = 2. From this we may use the recursion method to con- 


7 The new null region obtained in this way may be denoted by R(1, 2, 3, ; m.r ) and has the same size 


as R(1, 2. 3, ; m.c.) in figure 5. 


8See Pearson and Hartleys’ Tables (1941). 
® Evaluated by using inverse interpolation (see Whittaker and Robinson, 1934) in K. Pearson’s tables 
(1934a and 1934b). 
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struct the three-variate test as given in section 4.3, then from this, the 
four-variate test, and so on. 

The procedure of weighting for crowding and compensation for over- 
crowding will be seen to generalize fairly readily from the case n = 3 
to the cases n > 3. From this it seems reasonable to conclude that the 
average risk of the multiple comparisons test for the general case is fairly 
close to being a minimum. 


CONCLUDING REMARKS 


No attempt has been made in this paper to illustrate the various issues 
involved with examples from applied problems, or, at the end, to present 
the proposed test in a form convenient for practical applications. This 
has been done already in the preceding paper by the author (1951) and 
the reader is asked to refer to this for a discussion stressing the more 
applied aspects. 

The work in this and the preceding paper has been done with the 
financial assistance of a University of Sydney Pawlett Scholarship at 
the Iowa State College, a Commonwealth Research Grant at the Uni- 
versity of Sydney and a contract with the Bureau of Agricultural Eco- 
nomics at the Agricultural Experiment Station of the Virginia 
Polytechnic Institute. 
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Observations on the Egg Masses of Spotted 
Salamanders, Ambystoma maculatum 
(Shaw), in the Williamsburg Area 


Joun TuHornton Woon’ 
University of Virginia 


Rozert Houpen, WILKINSON 
College of William: and Mary 


Reported here are observations on the nesting of Spotted Sala- 
manders, Ambystoma maculatum (Shaw), on the Virginia coastal plain. 
The data on the time of egg deposition, the number of eggs voided in 
discrete groups, and the dimensions of egg masses are compared with 
published records of the nesting of this species in other regions. 

On February 13, 1949 nests of Ambystoma maculatum were dis- 
covered for the first time in Tidewater Virginia by Dr. J. T. Baldwin 
and Dr. B. M. Speese of the College of William and Mary. They 
collected several egg masses and adult specimens from roadside ditches 
bordering State Highway 5 about 2 miles west of Five Forks, James 
City County. Subsequent field work by the senior author disclosed a 
number of other breeding habitats in the Williamsburg area. One of 
these was noteworthy for its small water area and large number of egg 
masses. Excavation of a clay bank adjacent to the Confederate Redoubt 
one mile southeast of Williamsburg resulted in the formation of two 
pools. These pools are shallow (maximum depth, 12 inches), and have a 
combined surface area of less than. 300 square feet, yet in 1951 they 
contained 47 discrete egg masses. 

In any locality the date on which egg-deposition is initiated varies 
from year to year. Brimley (1939) reports that in the vicinity of Raleigh, 
North Carolina, breeding begins in late January in some years, and as 
late as March in others. This great latitude in the date on which egg- 
laying is initiated has not been confirmed by the observations of Wright 
(1908) and Bishop (1941), who found that in western New York, over 
a period of many years, the first egg masses of Ambystoma maculatum 
have invariably been noted between late March and mid-April. In the 
Williamsburg area eggs were first noted on February 13th in 1949, and 
on February 24th in 1951. 

In habitats near Williamsburg the egg-laying period of Ambystoma 
maculatum coincided with that of Hemidactylium scutatum in 1951, 
lasting from February 24th to March 10th. In some habitats both of 


1 Work on this report was completed when the senior author was a staff 
member of the Virginia Fisheries Laboratory and the Department of Biology 
of the College of William and Mary. 
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these species were encountered, but this was not usual. In a pond at 
Five Forks Hemidactylium females and eggs were found in the moss, 
Thuidium delicatulum (Hedw.) Mitt. along the banks, and masses of 
Ambystoma maculatum eggs were noted nearby attached to submerged 
sticks under from four to ten inches of water. 


The eggs of Ambystoma maculatum are deposited in a matrix of 
jelly. The.dimensions of an egg mass are partly a function of the number 
of eggs involved, but are also influenced by the amount of water which 
has been absorbed by the eggs and jelly, and by the extent of cohesion 
of the jelly matrix. Bishop (1941) noted that in New York egg masses 
in early stages of development had diameters of from 214 to 314 inches, 
but he pointed out that continued water absorption and partial disintegra- 
tion of the matrix may greatly increase the diameter, “. . . one such 
distended mass measured 414 inches in diameter.” Many Virginia egg 
masses were larger than those reported by Bishop. The smaller masses 
tended to remain spherical when suspended in water, but the larger 
aggregations of eggs appeared as elongated oval masses. The largest 
Virginia egg mass measured 156 mm. long, 119 mm. wide, and 43 mm. 
deep; the smallest mass measured 76 mm., 63 mm., and 28 mm. for the 
respective dimensions. 


Bishop (1941) found that the ovaries of a gravid Ambystoma 
maculatum may contain over 200 eggs, and that most nests contain less 
eggs than the probeble full overian complement of a female. Although 
some nests contained 200 eggs, others contained as few as 12, with the 
number generally nearer 100 (Bishop, 1941; Pope, 1944). Bishop (1943) 
reports that a female may deposit her eggs in one or two large masses, 
or in a half dozen small ones. These reports indicate that an egg-group 
found in the field is not likely to be correlated with the number of ovarian 
eggs in the gravid female that deposited it. The possibility that more 
than one female may deposit their eggs in a single nest is suggested by 
Hay’s (1893) report of egg masses containing from 300 to 400 eggs. It 
is also possible that this number of eggs was arrived at by estimate rather 
than through a direct count, or that aberrant specimens were involved. 


Twenty-one nests from the Williamsburg area were dissected to 
determine the numbers of eggs involved. These nests contained the fol- 
lowing numbers of eggs: 190, 164, 151, 151, 150, 129, 125, 125, 124, 122, 
113, 96, 89, 88, 81, 81, 77, 73, 57, 52. and 49. The mean of these 
counts is 109 eggs, with a standard deviation of 38.2. In nine New York 
nests Bishop (1941) found the following numbers of eggs: 256, 167, 
159, 144, 124, 120, 101, 75, and 61. The mean number of eggs in these 
nine nests is 134, and the standard deviation is 54.7. These data show 
that the number of eggs deposited by Ambystoma maculatum in a dis- 
crete mass is highly variable. Smith’s (1950) report that in Kansas 
each nest contains “approximately 100 eggs” is at variance with these 
findings. 
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SUMMARY 


In the Williamsburg area egg masses of Ambystoma maculatum 
were first noted in middle and late February, in 1949 and 1951. In the 
latter year, in the ponds near Confederate Redoubt located one mile 
southeast of Williamsburg, the last egg masses had been deposited by 
March 10th, the term of egg-laying lasting two weeks. Eggs were 
counted in 21 nests, and their numbers per mass ranged from 49 to 190. 
Little central tendency was evident from these data. Egg masses were 
variable in size and shape, with smaller ones usually appearing spheroi- 
dal, and larger ones frequently elongate and flattened. 
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HYPERKERATOSIS IN CATTLE 





The Production of Hyperkeratosis by the 


Administration of a Lubricant 
Witson B. Betu 


Virginia Polytechnic Institute 


Since the report on bovine hyperkeratosis (X-Disease) by Olafson 
(1947), observations of various workers have indicated that the con- 
dition is widespread in the United States. The disease has appeared on 
a number of farms in Virginia, and in some instances has caused con- 
siderable losses. The etiology of many outbreaks has remained obscure 
in spite of attempts to identify an infectious agent, a nutritional dis- 
turbance, or a toxic substance. During the investigation of an outbreak 
in March 1951 it became apparent that a toxic substance or substances 
might be the causative agent. Materials were obtained and experiments 
conducted to test this possibility. This report describes the results 
obtained with one of the materials, a lubricant, with which symptoms 
and lesions of bovine hyperkeratosis were experimentally produced in 
calves. 


LITERATURE 


The early attempts to demonstrate the etiological agent, recently 
reviewed by Olson and Cook (1951), either failed or resulted in the 
production of one of two conditions, both unlike the spontaneously 
occurring disease. The early attempts at experimental production of 
bovine hyperkeratosis failed to produce the characteristic lacrimation, 
skin lesions, or the typical chronic symptoms. Olson, Cook and Brouse 
(1950) reported on an outbreak in a group of calves fed various rations. 
They obtained data that indicated an unknown factor or factors in the 
particular alfalfa hay and dicalcium phosphate used that enhanced the 
severity of the disease. Olson and Cook (1951) later reported the develop- 
ment of hyperkeratosis in experimental calves fed the’ feed-stuff left 
over from the original feeding trial. The environment and contact were 
eliminated as possible factors. Olafson and McEntee (1951) were able 
to produce experimentally bovine hyperkeratosis by the feeding of a 
processed concentrate feed-stuff. Later McEntee, Hansel, and Olafson 
(1951) succeeded in producing hyperkeratosis in cattle with the free 
fatty acid fraction of the ether extract of the processed concentimte. 
Wagener (1951) published recently the results of his observations 
made from 1946 to 1948 on bovine hyperkeratosis in Germany. He was 
able to produce the disease with a particular wood preservative which 


Received for publication November 30, 1951. 
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affected the animals by contact, feeding, or inhalation. These recent 
developments indicate that bovine hyperkeratosis is caused by a toxic 
substance or substances, although the specific agent has not yet been 
reported. 


MATERIALS AND METHODS 


Approximately five pounds of lubricant (Lot 1) were obtained from 

a farm on which hyperkeratosis had occurred. Normal calves, ranging 

in age from 33 to 65 days, were obtained from the Virginia Polytechnic 

Institute herds and placed in a barn in separate stalls. No cattle had 

been housed in this barn for the past three years. The ration consisted 

of hay, grain, and milk obtained from the parent herd in which hyperkera- 
tosis is known not to have occurred. The calves were all handled in the 
same manner, except for exposure to the lubricant, and no attempt at 
isolation was made. Four calves were administered varying amounts of 
lubricant by means of gelatin capsules; two calves were given varying 
amounts of a lubricant (Lot 2) purchased on the market, and two calves 
served as controls. One control calf (No. 241) was kept in a separate 
stall while the other control (No. 1356) was kept for a total of 65 days 
in a stall with calves receiving the lubricant. The rate and duration of 
administration of the lubricant were varied according to the response of 
the calf. In all cases the administration was stopped when it became 
obvious that irreparable damage had occurred. 

The lubricant was administered according to the following schedules: 

Calf 242—Lot 1 lubricant. 15 grams daily for the first 7 days, none 
for the next 5 days, then 10 grams daily for the next 12 days 
of administration. 

Calf 1351—Lot 1 lubricant. 15 grams daily for the first 13 days, 10 
grams on the 14th day, then 15 grams on the 15th and 16th 
days. 

Calf 247—Lot 1 lubricant. 15 grams daily for the first 13 days, 10 
grams daily for the next 4 days, then 5 grams on the 18th 
day of administration. 

Calf 252—Lot 1 lubricant. 5 grams daily for the first 10 days, then 
5 grams on the 12th, 14th, 16th and 17th days. 

Calf 1366—Lot 2 lubricant. 10 grams daily for the first 12 days, none 
on the 13th day, then 10 grams daily for the next 5 days. 

Calf 1371—Lot 2 lubricant. 15 grams daily for the first 12 days, none 
on the 13th day, then 15 grams daily for the next 5 days. 


RESULTS AND DISCUSSION 
The four calves exposed by means of oral administration to the 


lubricant obtained from the farm and the two calves similarly exposed 
to the purchased lubricant developed symptoms and lesions associated 
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with bovine hyperkeratosis; no symptoms or lesions developed in the con- 
trols. During the period of study, July 10, 1951 to October 25, 1951 
and subsequently, no symptoms or lesions of hyperkeratosis and no 
evidence of a dietary deficiency appeared in 35 calves in the parent 
herd which received the same ration as the calves used in the experiment. 


Since the response of the calves exposed to the lubricant obtained 
from the farm was similar to the response of the calves to the purchased 
lubricant, the two groups will be discussed together. The first symptom 
noticed in all calves was lacrimation, which began from the 5th to the 
13th day and generally lasted throughout the duration of the disease. 
In some cases the lacrimation was so profuse that the cheeks remained 
wet. A slight rise in temperature (103.0°F. to 104.6°F.) lasting for a 
few days, occurred in four calves about the time of the appearance of 
lacrimation. All calves showed slight loss of appetite during the early 
phases of symptoms. A transient diarrhea appeared in three calves from 
the 5th to the 11th day, and in two calves in the later stages, on the 
27th and 29th days. Four calves showed some salivation, which began 
as early as the 9th day and as late as the 18th day. 


Red areas on the lips and hard palate occurred in all calves. These 
areas appeared on the 8th day after beginning the feeding in one case 
and as late as the 22nd day in another. These areas disappeared in some 
instances, to be followed by development of similar areas at other sites. 
The gums (Fig. 1) of four calves became red and swollen. This change 
began at one incisor and in some cases progressed until the entire gum 
was involved. Papilloma-like proliferative changes developed in each 
calf either on the tongue, gums, dental pad, or muzzle (Figs. 2, 3, 4, 5). 
One calf developed such changes at all of these sites and a second calf 
developed them at three of the areas. 


The skin of all calves exposed to the lubricant became dry, hard 
and thickened. The folds of the skin became especially prominent over 
the face, neck and shoulders. These changes began as early as the 11th 
day from the initial feeding in one calf and as late as the 29th day in 
another (Figs. 5, 7). 


The same general picture was presented by all calves on autopsy. 
In addition to the skin and proliferative changes already described, gross 
changes were found in the digestive and urinary tracts of the calves. 
Flat ulcers were present in the abomasum of 5 calves; in some instances 
the mucosa of the intestinal tract was thickened. The kidneys appeared 
enlarged and pale, and pale streaks were present in the cortex. Each 
liver presented light colored areas, varying in size, on the surface; one 
liver had an extra-hepatic cyst. These areas extended into the par- 
enchyma. The gall bladder was distended in most calves, and the wall 
appeared thickened in all cases. Small cyst-like nodules on the lining of 
the gall bladder were seen in 5 calves, 
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Examination of tissues obtained at autopsy revealed histological 
changes that are in agreement with the changes reported by the authors 
previously cited. An accumulation of keratinized material and keratini- 
zation of the hair follicles developed (Fig. 8) although these changes 
were not so extensive as those seen in more chronic cases. ‘Cystic dilata- 
tions of the collecting tubules in the kidney cortex with accompanying 
fibrosis were present in all calves (Tig. 9). Areas of degenerating liver 
cells (Fig. 10) with occasional cellular infiltrations, were seen in some 
liver sections. Some areas of fibrosis were also present, as were prolifer- 
ations of the bile ducts. Cystic dilatations of the glands of the gall bladder 
were observed in some instances (Fig. 11). These occurred at the site 
of the nodules observed at autopsy. The kidney changes were uniform 
for all calves; however, the other changes described do not represent a 
complete and detailed study of all individuals. 


The cause of the condition produced by the lubricant was rapid 
when compared to spontaneous cases of hyperkeratosis. One calf, No. 
1351, developed a secondary abscess of the tongue and pharynx and 
died on the 25th day. The remainder of the calves were allowed to go 
until they became extremely weak, at which time autopsies were per- 
formed. The duration, under these conditions, ranged from 27 to 39 
days. The amount of lubricant administered varied from 70 grams ad- 
ministered in 17 days, to 255 grams administered in 18 days. There was 
no recognizable relation of severity of the symptoms produced to the 
amount administered. The initial symptoms of lacrimation appeared in 
Calf 252 in 7 days and after a total of 30 grams of lubricant while 
in Calf 1351 lacrimation appeared on the 11th day after a total of 150 
grams of lubricant. There appeared no relation between the gross 
changes produced and the degree of exposure, except possibly in Calf 
252. This calf, which received the smallest amount, developed the most 
typical proliferative changes of the oral mucosa. This suggests that the 
administration of smaller amounts of the lubricant might result in a more 
prolonged and chronic syndrome. The cessation of administration of the 
lubricant resulted in no improvement in the condition of the calves, and 
it is apparent that the amounts given were more than sufficient to bring 
about the symptoms and lesions observed. The time of appearance of the 
first symptoms, subsequent developments, and the duration of the dis- 
ease produced are in general agreement with the data reported by 
McEntee, Hansel, and Olafson (1951). 


LEGEND FOR FIGURES 


Fig. 1—Gums reddened and swollen, Calf 247, at 21 days. Fig. 2— 
Oral papillomas, Calf 252, at 25 days. Fig. 3—Papilloma on dental pad, 


Calf 252, at 25 days. Fig. 4—Wart-like proliferation of muzzle, Calf 


252, at 25 days. Fig. 5—Section of proliferative lesion mouth, Calf. 252. 
Fig. 6—Hyperkeratosis, Calf 242, at 36 days. 
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SUMMARY 


Hyperkeratosis has been experimentally produced in calves by the 
administration of a lubricant recovered from a farm on which the disease 
had occurred and by a similar product purchased on the market. Typical 
symptoms and lesions were produced in six calves while two controls, 
one of which was in contact exposure, remained normal. As little as 70 
grams administered over a period of 17 days produced the symptoms and 
lesions. As far as the author is aware this is the first experimental pro- 
duction of bovine hyperkeratosis with a petroleum product. 
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Fig. 7—Skin area of neck, Calf 1371, at 38 days. Fig. 8—Section of 
skin, Calf 1366, showing keratinized layer and keratinization of hair 
follicles. Fig. 9—Section of kidney, Calf 242, showing cystic dilatation of 
tubules. Fig. 10—Section of liver, Calf 242, showing degeneration of 
—. Fig. 11—Section of gall bladder, Calf 1366, showing dilatation of 
glands. 
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News and Notes 


(Editor’s note: News contributions should be sent to the 
person whose name appears at the end of the appropriate 
section.) 


NOTES FROM THE PRESIDENT 


I have a few brief announcements and comments for which I beg a 
moment of your attention. 


Your Council, at the request of the Chamberlin Hotel at Old Point 
Comfort, has agreed to postpone our annual meeting to the third week- 
end in May. 1952, that is, the 15th, 16th, and 17th. This was done 
because the original date agreed upon proved difficult for the hotel 
management We hope this wili not inconvenience you and that you will 
make your plans accordingly. 

May I cell attention to the notices which you will receive from your 
Section Secretaries in February calling for papers for the May meetings. 
Please cooperate with the deadline which he sets for receiving these titles 
for he must make up your Section program and submit it to our Editor 
by March 5th. At the same time your Section Secretary will request you 
to submit an abstract of your paper not to exceed 200 words. Please 
give your attention to this matter because if it is too long, the Editorial 
Staff will find it impractical to consult with you individually and will 
be forced to resort to their own judgment with regard to its reduction. 

At the suggestion of Mr. John Thornton Wood, Chairman of our 
Committee or Resource-Use Education, I wrote Governor John S. Battle 
requesting him to examine the possibility of making available to the pub- 
lic schools literature on Conservation in a form that would be of value 
to these groups. He referred this matter to the Department of Education 
and it is our hope that steps in this direction will be taken. 

The Administrative Secretary of the American Association for the 
Advancement of Science requested the members of the Association’s 
Council to write to the Senate Appropriations Committee to restore the 
‘budget of the National Science Foundation reduced by the House Com- 
mittee from the requested fourteen million dollars to three hundred thou- 
sand dollars. Nearly thirteen millions of the requested budget had been 
earmarked for basic research and for a fellowship program. I wrote in 
the name of our Acedemy to the Chairman of the Senate Appropriations 
‘Committee. I am recently informed that budget was finally approved at 
three million dollars. I regret that the letters of scientists did not move 
the Congressmen to restore more of the requested budget. 

Dr. Henry Leidheiser, Jr. has written a fine article in which he 
tries to interest more industrial scientists in our State in the Virginia 
Academy of Science and in presenting more papers at our annual meet- 
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ings. (See The Bulletin of the Virginia Section of the American Chemical 
Society, Vol. 29, numbers 1 and 3, 1951.) I wish to extend a cordial 
invitation to all of these scientists not yet members to join this Academy 
and take an active part on its programs, 

Finally. I would request you to send me such suggestions as you 
deem appropriate and feasible of execution designed to make our Academy 
function more effidiently for the cause of scientific research, science 
teaching, and for the development of interest in science in Virginia. 

Paut M. Parrerson 


RESEARCH COMMITTEE 


Encouragement of research in Virginia has always been one of the 
chief objectives of the Academy. Research Grants-in-Aid, and also the 
annually available J. Shelton Horsley Research Award, are offered by 
the Academy to encourage and stimulate research work of its members. 
The attention of Academy members is called to the availability of both 
the Grants-in-Aid, znd the Award. The Research Committee urges the 
greater interest and participation of members in both. 


ResEARCH GRANTS-IN-AID 


The income from the Academy’s Endowment Fund is set aside each 
year to provide assistance to Academy members who submit to the 
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Research Committee worthy research problems for which financial aid 
is needed. Research grants have varied in amount but during recent years 
have averagec approximately $100.00 each. 


The appiicant is requested to outline his problem and list the items 
for which assistance is required. The names and addresses of three 
persons who know his qualifications for scientific research should be 
submitted with the request. 


Grants will be made on the basis of (1) apparent merit of the 
research project, and (2) the need of the grant to aid in carrying out 
a project, or some portion of it. Worthy problems being worked on by 
younger members, and by members on staffs of smaller schools—as well 
as investigations by older members and by those on staffs of larger 
institutions—will be given full consideration in the awarding of grants. 


Since funds are now available, the committee urges that applications 
be mailed at an early date to W. S. Flory, Jr., Boyce, Virginia, so that 
they may be -onsidered at its next meeting. 


J. Sue.ttron Horstey Researcu AwWArpD 


Members of the Academy and officers of Sections are requested to 
give thought to the entering of papers for the 1952 Horsley Research 
Award. The only requirements for eligibility for the J. Shelton Horsley 
Research Award in 1951 were that submitted papers should (1) present 
original research (2) by Academy members (3) at the annual May 
meeting, - (4) the paper dealing either with unpublished data or be a 


published work - if publication had not been prior to May of the previous 
year. 


It is anticipated that essentially the same rules will apply to eligi- 
bility in 1952 as in 1951. Details of these will be worked out by the 
Research Committee early in 1952, and will be supplied to every member 
by the secretaries of the varicus Sections with their calls for papers 
for the May meeting. Members desiring details of eligibility before 
receiving the secretaries call for papers will be furnished with this 
information or inquiry to the Research Committee. 

The Section Officers are again asked to serve as advisors to the 
Research Committee and to pass on to it papers which they deem prize- 
worthy, even though they are not formally submitted for the award. This 
procedure in 1950 and 1951 resulted in a marked increase in the number 
of papers considered. 


The Horsley Award is one of the top prizes available to scientists 
in Virginia. 'The widespread interest and participation of members in 
the seeking of this award will mark a healthy condition with respect to 
the type of research work being carried out in the state. Such interest and 
participation is desired and solicited ——W. S. Fiory, Jr., THe Buanpy 
EXPERIMENTAL Farm, University or Viretnia, Boyce 
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INSTITUTE OF TEXTILE TECHNOLOGY 


Eleven Research Fellows arrived at the Institute of Textile Tech- 
nology, Charlottesville, on September 13 to begin their first year in the 
Institute’s two-year graduate program. 

Their arrival marks the fifth academic year of the Institute edu- 
cational program which leads to che Master of Science degree in textile 
technology. Only students showing a special interest and aptitude in 
textile science or having practical experience in textile work are accepted 
by the Institute. 


The new students are: Richard P. Barber, Bradford Durfee Techni- 
cal Institute: Anibal L. Ferreira, New Bedford Textile Institute; Isaac 
Jarkowsky, Brooklyn College; Lawrence C. Legere. New Bedford Tex- 
tile Institute; Earl J. Nickerson, Lowell Textile Institute; Arthur J. 
Pendleton, Rhode Island Scheol of Design; Robert G. Seid, Georgia 
Institute of Technology; Robert E. Staples, Bradford Durfee Technical 
Institute; Albin F. Turbak, New Bedford Textile Institute; James W. 
Whitworth, Lowell Textile Institute; Theodore Williamson, Jr., Brad- 
ford Durfee Technical Institute. 


AGRICULTURAL SCIENCE SECTION 


Agricultural workers in Virginia regret to hear that Mr. V. A. 
Tiedjens has resigned as Director of the Virginia Truck Experiment 
Station to accept a position with a fertilizer concern in Ohio. Mr. Tiedjens 
has served the vegetable growers faithfully and well during his years at 
the Norfolk Station. Although the horticultural industry in Virginia is 
suffering a real loss in the resignation of Mr. Tiedjens all his friends 
join in wishing him the best of saccess in his new work in Ohio. 

The new Director at the Virginia Truck Experiment Station will be 
Mr. William H. Brittingham who has been horticulturist at the Texas 
Agricultural and Mechanic Arts College since 1944. Mr. Brittingham is 
a plant breeder by training. He was born in Portsmouth, Virginia in 
1907, secured his B.S. from Penn State and an A.M. degree from Prince- 
ton in 1929 His Ph. D. degree was obtained at Maryland in 1942, with 
a major in plant breeding. 

Mr. Flood Andrews has been granted a leave of absence from his 
duties in thc vegetable section of the department of horticulture at 
Virginia Polytechnic Institute to assist in the Foreign Agricultural Pro- 
gram of the U.S. Department of Agriculture. Mr. Andrews will be 
stationed in Bolivia, South America. He will assist in setting up research 
and demonstration projects in vegetable production—W. P. Jupx1ns, 
Vireinta AcricutturaL Experiment Station, Biackssure. 


Astronomy, Mai1nematics, Puysics SEcTIon 


A new building to house the departments of Mathematics and 
Physics is under construction at Randolph-Macon College. It is to be 
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completed at the end of the present session. 


Mr. E. E. Gardner. former instructor in Physics at Virginia Poly- 
technic Institute has accepted a position as Assistant Professor of Physics 
at the Naval Academy. 

Mr. John Linsley has joined the staff of Rouss Physical Laboratory 
as Assistant Professor of Physics. He has recently been associated with 
Professor E P. Ney of the University of Minnesota, former member 
of this section, in cosmic ray research. 

Mr. A. R. Kuhlthau has returned to the University of Virginia as 
a member of the staff of the Naval Ordnance Laboratory. He has been 
on the Physics staff at the University of New Hampshire. 


Mr. J. L. Young, III has resigned from the staff of the Naval Ord- 
nance Laboratory at the University of Virginia to accept a position at 
Los Alamos. 

Mr. R. E. Alley, Jr., formerly at University of Richmond, is now 
a member of the research depariment of Bell Laboratories. Professor J. J. 
Taylor has been named Acting Chairman of the Physics Department in 
his place. 

New members of the Physics staff at University of Richmond are 
G. P. Williams, Jr. and C. O. Alley, Jr. Mr. Williams was a member 
of the Carson-Newman faculty before coming to Richmond. Mr. Alley 
returns from Princeton University where he was doing graduate work. 

Professor B. V. English of the Physics Department of Randolph- 
Macon College has been granted a year’s leave of absence to pursue 
graduate studies at Pennsylvania State College. 


Professor J. W. Beams of the University of Virginia has been 
appointed a member of the committee of the National Heart Institute. 

Three new members of the Physics staff at Virginia Military Insti- 
tute are M. J. Clingan, D. R. Carpenter, and J. B. Breazeale. Mr. Clingan 
was last year a member of the faculty of Newberry College. Mr. Car- 
penter comes from Cornell University and Mr. Breazeale from the 
University of Alabama.—I. G. Foster, Virginia Military Institute. 


BactTEerioLocy SEcTION 


The fall meeting of the Bacteriology section was held at the College 
of William and Mary cn November 17, 1951. After a luncheon at the 
Williamsburg Lodge the following scientific program was presented. 

A Simple Medium for the Detection of Pseudomonas, J. M. Sharpley. 
Laboratories, Dept. of Health, City of Richmond, Virginia. 

The Isolation of Enterococci from Urine, A. L. Rosensweig and 
Norman Ikari. McGuire Veterans Hospital, Richmond, Virginia. 

Types of Streptococci found in the Vagina of Mares, Melahat Pusat 
and P. Arne Hansen. Dept. of Bacteriology and Live Stock Sanitary 
Service Laberatory, University of Maryland, College Park, Maryland. 
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A Possible Source of Error in the Microscopic Examination of 
Smears for Tuberculosis, A. L. Rosensweig. McGuire Veterans Hospital, 
Richmond, Virginia. 

Streptococci as Secondary Invaders in Blackhead Liver Lesions, 
Arthur P. Harrison and P. Arne Hansen. Dept. of Bacteriology and 
Live Stock Sanitary Service Laboratory, University of Maryland, Col- 
lege Park, Maryland. 

Fermentation of Inositol by Propionibacterium pentosaceum, 
Wesley A. Volk. Dept. of Microbiology, University of Virginia, School 
of Medicine. Charlottesville, Virginia. 

Factors Relative to Growth of Endamoeba hystolitica, E. Clifford 
Nelson. Dept. of Bacteriology and Parasitology, Medical College of 
Virginia, Richmond, Virginia. 

Infectious Hazards of Bacteriological Techniques, Film by U. S. 
Public Health Service. C.D.C. 


IN MEMORIAM 


With the death of Doctor George McLean Lawson on September 
20, 1951, the members of the Bacteriology section of the Virginia 
Academy of Science lost one of their most highly esteemed members. 
Doctor Lawson’s leadership was instrumental in the development of 
this section of the Academy and with his passing the members will deeply 
miss a fine friend and wise counselor. 


—J. Doveras Rem, Mepicat CoLuece or VirGinia. 


Biotocy SEctTIon 


Mr. Austin B. Clark’s manuscript entitled “The Butterflies of Vir- 
ginia” has gone to the printer and should appear before 1951. There will 
be 31 plates illustrating each of the 155 species known from Virginia, 
and few that should occur here but which have not as yet been found. 

A grant-in-aid by the Virginia Academy in 1941 has led to the 
recent appearance of the first genethnics manual by E. S. Craighill 
Handy and Elizabeth Green Handy. “The Genethnic Screen” is pub- 
lished by Genethnics Incorporated, Oakton, Virginia. 

Mr. Colgate W. Darden, Jr., President of the University of Virginia, 
has recently announced the promotion of Mr. Edwin M. Betts from 
Associate Professor to Professor of Biology and Mr. B. F. D. Runk 
from Assistant to Associate Professor of Biology. 

Mr. Dabney Lancaster, President of Longwood College, has an- 
nounced the promotion of Mr. Robert T. Brumfield from Associate Pro- 
fessor to Professor of Biology. 

The following biologists from Virginia worked during the summer 
at Oak Ridge National Laboratory: 

Mr. Robert T. Brumfield 
Mr. James Norman Dent 
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Mr. Orland E. White, Director of the Blandy Experimental Farm, 
Boyce, Virginia has returned to the University of Virginia after serving 
the past year as Visiting Professor at University College, Mandalay, 
Burma.—Lanitey Hustep, University oF VirGINIA 


CuHeEMIstRY SECTION 


Frank Berton Carpenter. who was for fifty years Chief Chemist of 
the Virginia-Carolina Chemical Corporation, died on October 3, 1951, 
at the age of 87. He was one of the organizers of the Virginia Chemists’ 
Club, in 1907, and was later one of the chairmen of the Virginia Section 
of the American Chemical Society. 

Mr. Earl K. Fischer, former chief of the paint section of the Na- 
tional Bureau of Standards, who for four years was head of the physical 
chemistry division of the Institute of Textile Technology, at Charlottes- 
ville, died of a heart attack at his home in Kensington, Maryland, in 
August, 1951. 

From the Medical College of Virginia: Mr. Charles C. Clayton. 
assistant professor of biochemistry, has been awarded a grant of $5,184.00 
by the United States Public Health Service to apply to research on the 
relation of minerals to cancer. An additional $4,500.00 is ear-marked for 
the continuation of this study in its second year. 


Mr. Paul S. Larson, research professor, of biochemistry, and also 
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pharmacologist for the American Tobacco Company Research Laboratory, 
presented a paper at the American Chemical Society Jubilee in New York 
in September on “Present Knowledge Concerning the Metabolism of 
Nicotine by the Animal Body and the Nature and Action of Tobacco 
Smoke Irritants.’”” Later in September he presented a related paper at 
the World Tobacco Congress in Amsterdam, Holland. 


Mr. R. Blackwell Smith, dean of the School of Pharmacy and asso- 
ciate professor of pharmacology, presented a paper entitled ‘“Toxicologi- 
cal Studies on the Zinc and Disodium Salts of Ethylene-bis-dithio- 
carbamates” at the XIIth International Congress of Pure and Applied 
Chemistry in New York City in September. Associated with Mr. Smith 
in this work were Messrs. J. K. Finnegan, P. F. Sahyoun, H. B. Haag, 
and Paul S. Larson. 


The July number of the Proceedings of the Society for Experimental 
Biology and Medicine contained a paper by Messrs. G. Watson James, 
III, and Lynn D. Abbott, Jr., on the “Effectiveness of Thiocyanate Ana- 
log of Vitamin Bz in Pernicious Anemia.” 

Mr. J. C. Forbes and Mr. G. M. Duncan are authors of a paper on 
“The Effect of Acute Alcohol Intoxication on the Adrenal Glands of 
Rats and Guinea Pigs,” appearing in the September issue of the Quarterly 
Journal of Studies on Alcohol. 

From the Virginia Institute for Scientific Research: Mr. Carl J. 
Likes, former head of the department of chemistry at Hampden-Sydney 
College, and Mr. Victor J. Kehrer, formerly of the Nickel Fellowship 
at the Mellon Institute, have joined the staff of the Institute. 

From the University of Richmond: Mr. Thomas C. Franklin (Ph. D., 
Ohio State University) is assistant professor of chemistry in charge of 
physical chemistry. Mr. Lemuel F. Smith, professor emeritus of Kalama- 
zoo College, is visiting lecturer in chemistry. Puryear Fellowships for 
graduate study in chemistry have been awarded to Mr. Fletcher B. Owen, 


from the University of Richmond, and Mr. Ray D. Sothern, from Ursinus 
College. 


From the University of Virginia: Mr. Fred W. Young has returned 
from Oak Ridge as research associate in chemistry. Mr. Lewis B. Johnson 
is research associate in chemistry at the School of Engineering. Mr. Tom 
Crowell has returned from du Pont to his position at the University. 
Mr. John Cathcart has joined the Division of Metallurgy at Oak Ridge. 
Mr. Carl D. Lunsford has been awarded a predoctoral du Pont Fellowship 
for work under the direction of Mr. Robert E. Lutz. 

From the Virginia Polytechnic Institute: Under the direction of 
Mr. Nelson F. Murphy, a listing of “Chemical and Chemical Engineering 
Journals Available in Virginia” has been completed. 

Mr. Frank A. Vingiello spent the summer studying moderators for 
the Aircraft Nuclear Reactor at the Oak Ridge National Laboratory. 
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Mr. Russell E. Leed has been granted leave of absence for a year for 
work with the Atomic Energy Commission at Oak Ridge. Col. W. L. 
McPherson participated in the 23rd Radioisotopes Techniques Course 
during the summer at Oak Ridge. 


From Hollins College: Miss Lois Ann Trafton (A. B., Hollins) of 
Norfolk has joined the staff of the department of chemistry. 
—Witu1a E. Trovt, Jr., Box 216, 
University or RicHMoND, VIRGINIA 


ENGINEERING SECTION 


The Central Virginia Engineers’ Club is an organization of people 
engaged in or interested in the profession and practice of Engineering, 
bound in fellowship by the desire for increased professional unity and 
social relaxation. 


Its membership is made up principally from residents of the Rich- 
mond, Petersburg, and Hopewell areas. Meetings are held jointly with 
one of the National Societies such as American Society of Mechanical 
Engineers, American Society of Civil Engineers, American Society of 
Refrigerating Engineers, American Institute of Electrical Engineers, 
American Institute of Architects, and American Chemical Society. In 
addition to the meetings several inspection trips have been made. The 
membership is approximately four hundred. 
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The Department of Chemical Engineering, Virginia Polytechnic 
Institute is carrying on research on the electrodeposition of aluminum, 
financed by a grant from the Air Material Command. 


Fred W. Bull, Associate Professor of Chemical Engineering, Vir- 
ginia Polytechnic Institute, has been on leave, for work at Oak Ridge. 

V. G. Szebehely, Applied Mechanics Department, Virginia Poly- 
technic Institute., resigned in September, 1951, for a research position 
with the David Taylor Model Basin. 

Dan Frederick returned to Virginia Polytechnic Institute after a 
year’s leave of absence at the University of Michigan, where he completed 
all requirements for his Ph. D. in Mechanics except his thesis. This he 
is now working on. He was promoted to an Associate Professor in Ap- 
plied Mechanics in September, 1951. 

D. H. Pletta, Head of Applied Mechanics, Virginia Polytechnic 
Institute, attended the annual meeting of the ASCE in New York, 
October 22-26, for the dedication ceremonies of the newly formed En- 
gineering Mechanics Division. He is Chairman of the Technical Com- 
mittee on Experimental Analysis and Analogues. 

Mr. Pletta also attended the annual convention of the Virginia 
Society of Professional Engineers at Old Point Comfort on November 
1-2. He is First Vice President of this society for the coming year, and 
is also a member of the Committee on Education of the National Society 
of Professional Engineers. 


The Virginia Society of Professional Engineers, Inc., Mr. John A. 
Rives, President, has started publication of a magazine of news and 
articles relating to professional engineering. Called The Virginia En- 
gineer, the journal is published quarterly under Editor H. H. Roberts. 


The Virginia Society of Professional Engineers held their third 
annual meeting at Old Point Comfort on Nov. 2-3, 1951. The program 
included talks by Louis Lee Guy, Member of the Virginia House of 
Delegates on “Is Legislation Desirable”; by T. Keith Legare, Executive 
Secretary of the National Council of State Boards of Engineering Ex- 
aminers, on “The Development of Engineering Registration”; and by 
L. L. Dresser, President of the National Society of Professional En- 
gineers, on “Professional Aspects of Members and the National Society.” 

Mr. Frank C. Vilbrandt, Head of the Department of Chemical 
Engineering at Virginia Polytechnic Institute, was an Executive Staff 
member of the Temporary National Office of Ordnance Research, De- 
partment of the Army, located at Duke University, Durham, North Caro- 
lina, from August 13 to September 15, and is being retained as consultant 
of the permanent organization, but he is now back in Blacksburg. 

Mr. N. F. Murphy, Research Professor in Chemical Engineering, 
was appointed as delegate to the World Metallurgical Congress in Detroit, 
Michigan, October 14 to 20, where he was a member of a group of 
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international scicntists assembled for the exchange of information. 
Mr. Murphy was assigned to the specific group which discussed the 
various problems in the field of metallurgical education. 

—Netson F. Murpny, Vireinia Potytecunic Institute 


PsycHotoey SEcTION 


A number of State psychologists attended the annual meetings of 
the American Psychological Association in Chicago, August 30 - Septem- 
ber 5. From the University of Virginia Mr. Frank A. Geldard spoke at 
a symposium: Relations between Governmental and Non-Governmental 
Development of Psychological Resources, and was chairman of one of 
the sessions sponsored by the Division of Military Psychology, of which 
he is president. Mr. Richard H. Henneman presented a paper, A Com- 
parison of Vision and Audition as Sensory Channels for Communication. 
Mr. L. Starling Reid read a paper, The Development of Non-Continuity 
Behavior through Continuity Learning. Several psychologists from the 
Lynchburg State Colony also attended the APA meetings: Mr. John N. 
Buck, Miss Hannah S. Davis, Mr. Allen Cohen, and Mr. William Dakos. 
Mr. Buck was chairman of a round table: Theoretical and Practical 
Aspects of the H-T-P Technique; Miss Davis presented a paper, The 
Case of the Rays. 


Mr. John N. Buck retired in September as chief psychologist of the 
Lynchburg State Colony. His present address is P. O. Box 157, White 
Stone, Lancaster County, Virginia. Miss Hannah S. Davis was appointed 
acting chief psychologist to succeed Mr. Buck at the Colony. Miss Lois 
Brinkmann, of Miami, Florida, was appointed as a psychological interne 
at the Colony on October 1. 

At a meeting of the State Board for the Examination of Clinical 
Psychologists held at the Lynchburg State Colony on August 4,. 
Mr. Reuben Horlick was certified as a clinical psychologist in the State 
of Virginia. 

The Advisory Psychophysiology Panel of the Office of Naval Re- 
search met at the University of Virginia on November 12 and 13. 
Mr. Frank A. Geldard, a member of the panel and director of a Naval 
research contract on tactual vibratory sensitivity at the University was 
host to the group. 


The Section Secretary has been notified by the Executive Secretary 
of the American Psychological Association of the elimination of the 
membership class of State Affiliate. This action by the Council of Repre- 
sentatives of the APA will in effect discontinue the status of membership. 
in the APA of several members of the Psychology Section of the Vir- 
ginia Academy of Science. The membership status of State Affiliates was 
terminated as of December 31, 1951, along with the attendant privilege 
of subscribing to APA journals at the special rates. Those individuals 
interested in becoming associate members of the APA may obtain infor- 
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mation relative to qualifications and procedures from either Mr. S. B. 
Williams of the College of William and Mary, or from the Executive 
Offices of the American Psychological Association, 1515 Massachusetts 
Avenue, N. W., Washington 5, D. C. 


A course in the Rorschach method is being planned at the University 
of Virginia through the Department of Education for the Spring Semes- 
ter. The course will be given in three terms, the first year being devoted 
to basic principles of scoring and interpretation, with the third term 
being concerned with advanced interpretation. It will be open to practic- 
ing clinicians and psychologists with advanced graduate training. For 
further information applicants may write to Lindley J. Stiles, Dean, De- 
partment of Education, University of Virginia, or to the Instructor, 
Arthur J. Bachrach, Director, Division of Clinical and Medical Psy- 
chology, University of Virginia Hospital, Charlottesville, Virginia. 

—Ricnarp H. Henneman, UNIversity oF VIRGINIA. 


Sratistics SECTION 


Mr. Charles Mottley, Mr. J. W. Youden, Mr. Boyd Harshbarger,. 
Mr. R. A. Bradley, Mr. G. W. Tyler, Mr. D. B. Duncan, Mr. M. E. 
Terry attended the joint Institute of Mathematical Statistics and 
American Society for Quality Control meetings in celebration of the 50th 
Anniversary of the National Bureau of Standards. Mr. D. B. Duncan 
presented a paper, and Mr. J. W. Youden an invited address. 

The Department of Statistics of Virginia Polytechnic Institute pre- 
sented as its fall quarter invited speakers: Mr. John W. Tukey, Prince- 
ton, November 2 and 3; Mr. John H. Nair, Director of Research of 
Thomas Lipton Company, November 26; Mr. D. B. DeLury, Ontario 
Research Council, December 7 and 8. The department has fourteen gradu- 
ate students with Mr. Thomas Russell, teaching fellow and Mr. Clyde 
Kramer in charge of the Statistical Laboratory. Two semi-annual reports 
on research on statistical methods for sensory difference testing are now 
available from the Department of Statistics of Virginia Polytechnic 
Institute. 

Mr. Charles Beasley, Fellow of the Institutes of Public Health has 
resigned his fellowship on the activation of his reserve commission, United 
States Army. 

Mr. E. L. Cox is now on the staff at Dugway Proving Ground, Utah. 

The Institute of Mathematical Statistics and the E.N.A.R. section 
of the Biometrics Society will hold a joint Spring meeting March 19-21, 
1952, on the campus of Virginia Polytechnic Institute. Dormitory facili- 
ties will be provided for male graduate students for $1.00 per night. 

Mr. Elie Weeks and Mr. John W. Griswold, Quartermaster Board, 
Fort Lee visited the Department of Statistics on Oct. 9. Mr. Griswold 
has joined the Academy. 
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Members of the section are invited and encouraged to present papers 
at the May meeting of the Academy. Titles should be sent to the secretary 
no later than February 18, 1952.—Murton E. Terry, Virarnia Poty- 
TECHNIC INsTITUTE. 








